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SUPPRESSION  OF  ICE  FOG  FROM  COOLING  PONDS 


by  .j 

Terry  McFadden 

INTRODUCTION 

Icc  fog  is  caused  by  extremely  small  ice  crystals  (8*/rr  to  3 Sum  diameter)  dispersed  in  the  air. 

It  forms  whenever  temperatures  fall  below  approximately  -30aC  and  water  vapor  is  introduced 

into  the  air  in  excessive  quantities.  ! 

Several  studies  have  delved  into  the  nature  and  causes  ot  ice  fog.  Benson  (1965)  made  the  first 
comprehensive  study  into  the  makeup  of  the  phenomenon  and  estimated  that  power  plant  cooling 
water  produces  over  60%  of  the  total  ice  fog  in  the  Fairbanks.  Alaska,  urban  area.  Ohtake  (1970) 
studied  the  nature  of  the  fog  crystals  and  the  mechanisms  of  fog  formation.  Suppression  of  ice  fog 
was  largely  ignored  until  McFadden  (unpublished  work)  designed  and  tested  a small  exhaust-gas 
dehydration  unit  in  1968.  In  1969  Tedrow*  (personal  communication)  designed  and  tested  a suc- 
cessful exhaust  dehydrator  on  an  Army  2'  -ton  truck.  Furnace  stack  dehydrators  had  been  used 
sporadically  with  mixed  results,  until  finally  McKay  ♦ (personal  communication)  designed  and  tested 
a stack  dryer  for  a boiler  in  1971 . Suppression  of  emissions  from  power  plant  cooling  water  had  not 
been  approached  prior  to  this  study. 

In  1969  the  Alaskan  Air  Command,  recognizing  the  problems  associated  with  ice  fog  at  Eiclson 
Air  Force  Base,  asked  the  Arctic  Health  Research  Center  to  investigate  ice  fog  suppression  on  the 

base.  It  was  proposed  that  the  major  contributor  to  fogging  problems  be  considered  first,  and  j 

- knowledge  gained  in  this  effort  could  then  be  used  to  formulate  attacks  on  the  remaining  sources.  j 

The  purpose  of  this  study  was  to  test  ideas  for  ice  tog  suppression  on  the  Eiclson  AFB  power  plant 
cooling  pond  and  to  report  on  the  results. 

The  Eielson  AFB  power  plant  cooling  pond  is  located  on  the  east  side  of  the  base.  Terrain  to 
the  east  of  the  pond  rises  to  an  altitude  of  500  m MSL.  The  pond  is  168  m wide  and  305  m long, 
and  is  divided  into  two  units,  76  m wide  and  305  m long,  by  a 1 5-m-wide  dike.  Surface  area  is 
46,360  mJ.  Depth  is  3 m and  very  uniform. 

Power  plant  cooling  water  at  approximately  25°C  is  injected  into  the  north  pond  (which  has 
been  designated  as  the  receiving  pond)  and  drawn  out  of  the  south  (intake)  pond  at  1 5°C.  Flow 
through  the  pond  is  uniform  at  0.38  m3  s'1 . Heat  input  to  the  pond  averages  1 8x  107  W. 

During  the  winter  of  1969  a preliminary  study  was  conducted  to  gain  background  information 
on  the  pond  and  its  operation.  A dewatering  pump  was  used  to  divert  the  water  from  the  receiving 
valve  house  to  the  intake  side  of  the  pond.  This  was  found  to  be  a feasible  means  of  isolating  the 
receiving  side  of  the  pond  so  that  it  could  cool  and  form  an  icc  cover.  Short-circuiting  the  pond 
in  this  manner  did  not  appear  to  pose  any  significant  overheating  problem,  and  plans  were  made 
to  conduct  a full-scale  experiment  during  the  following  winter  season. 

■i 

i 

7 *J.  Tedrow,  Mechanical  Supervisor,  USA  CRREL  Alaskan  Projects  Office.  j 

I fR.  McKay,  Professor  of  Mechanical  Engineering,  University  of  Alaska.  j 
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Figure  1.  E/e/son  A FB  power  plant  cooling  pond. 


The  winter  of  1970-71  was  one  of  the  mildest  on  record.  Temperatures  did  not  drop  sufficiently 
to  permit  any  worthwhile  experimentation,  and  work  on  the  project  was  postponed  until  the  winter 
of  1971-72.  Plans  were  made  to  install  a flood  control  pump  to  divert  the  water  from  point  1 (see 
Fig.  1 ) to  point  2,  thus  short-circuiting  the  normal  flow  pattern  and  isolating  the  receiving  side  of 
the  pond  completely. 

In  the  summer  of  1972  the  required  installations  were  completed  and  the  experiment  was  readied 
for  the  coming  winter  season.  During  the  winter  of  1972-73,  the  project  was  operated  on  a 24-hour 
per  day,  seven  days  per  week  basis,  from  early  December  to  mid-March.  The  entire  pond  was  frozen- 
over  and  ice  fog  suppression  was  achieved,  while  melting  ice  consumed  the  excess  heat  load  of  the 
power  plant.  A heat  balance  was  determined  for  the  pond,  and  radiative  and  evaporative  heat  losses 
were  measured.  Various  spray  and  flooding  techniques  were  investigated  for  their  efficiency  in  ice 
building,  and  freezing  rates  were  measured.  Evaporation  suppression  was  measured  while  using 
monomolecular  films  of  hexadecanol,  and  also  while  using  ice  cover.  Approximately  3.6x  10*  kg 
of  ice  was  formed  on  the  pond  and  subsequently  melted  by  the  cooling  water  from  the  power 
plant.  Heat  transfer  from  pond  water  to  the  ice  sheet  was  monitored,  and  heat  transfer  coefficients 
were  determined  for  comparison  with  calculated  values.  Temperature  and  wind  velocity  profiles 
above  both  open  water  and  ice  were  measured  to  establish  the  film  coefficients  of  heat  transfer  in 
these  cases. 

Section  I of  this  report  is  a short  review  of  ice  fog,  its  causes  and  effects.  Sections  II,  III  and  IV 
present  a discussion  of  the  various  components  of  the  heat  budget  equation.  Section  V reviews 
suppression  techniques  and  experiments  in  suppression  using  monomolecular  films  and  Section  VI 
presents  the  experiments  in  ice  fog  suppression  using  ice  cover.  In  Section  VII  recommendations 
and  conclusions  are  offered. 
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SECTION  I.  ICE  FOG 


Ice  fog  is  a phenomenon  peculiar  to  locations  where  temperatures  drop  below  -30°C  and  water  j 

sources  are  available.  (Water  sources  such  as  hot  springs,  open  water  on  rivers  or  lakes,  products  of  .V 

combustion,  and  even  breathing,  all  may  contribute  to  ice  fog.)  Ice  fog  differs  from  water  fog  in  < 

that  it  is  composed  of  very  small  ice  crystals  (8  jum  to  35  ion  diameter)  rather  than  water  droplets. 

Ice  crystals  of  this  size  are  very  slow  to  precipitate  (precipitation  rates  of  5x  10"4  to  5x  10"6  g HjO  ■ 

m~J  s'1  are  reported  by  Benson  1965),  and  thus  they  remain  in  the  air  for  long  periods  of  time.  \ 

The  size  of  the  particles  is  a function  of  the  cooling  rate  of  the  water  droplets  from  which  they  are 
formed. 

The  major  contributors  to  ice  fog  in  the  Fairbanks  area  have  been  identified  by  Benson  (1965) 
as  open  water  from  power  plant  cooling  requirements,  combustion  by-products  from  power  plants 

and  homes,  and  exhaust  emissions  from  automobiles.  Of  these,  he  estimates  that  open  water  from  ■ 

the  power  plant  cooling  ponds  of  Fort  Wainwright,  and  open  water  on  the  Chena  River  created  by 

the  Fairbanks  Municipal  Utilities  System  cooling  water  discharge  are  responsible  for  over  half  the  1 

ice  fog  in  the  Fairbanks  area.  Eielson  Air  Force  Base,  located  22  road-miles  (35.4  km)  southeast  j 

of  Fairbanks,  is  analogous  to  the  Fairbanks-Ft.  Wainwright  area,  but  of  somewhat  smaller  scale.  It  j 

differs  in  a few  major  respects,  however,  and  these  should  be  discussed.  j 

The  power  plant  stacks  at  Eielson  Air  Force  Base  are  37  m high  and  they  discharge  their  gases  at  t 

a relatively  high  velocity.  This  generally  has  the  effect  of  putting  these  gases  above  the  critical  level, 
so  that  they  do  not  merge  with  the  ice  fog  blanket  on  the  ground  within  the  confines  of  the  Air 
Force  Base.  They  are  not  normally  a major  contributor  to  ice  fog  on  the  base  (Robinson  1963). 

Only  under  certain  conditions  of  prolonged  ice  fog,  when  the  temperature  inversion  layer  has 
thickened,  will  gases  from  the  power  plant  stacks  be  forced  back  into  the  ice  fog  blanket.  This 
situation  rarely  occurs  within  the  base  boundary.  However,  meteorological  conditions  can  cause 
the  suck  gases  to  be  a major  contributor  to  ice  fog  in  the  surrounding  residential  areas. 

A large  percentage  of  the  heating  at  Eielson  AFB  is  done  with  steam  from  the  power  plant  (in 
contrast  with  Fairbanks  where  condenser  heat  is  largely  discharged  into  the  Chena  River).  This  has 
the  effect  of  decreasing  the  ice  fog  in  the  housing  areas,  since  buildings  heated  by  power  plant 
steam  are  not  emitting  products  of  combustion  from  individual  furnaces.  The  cooling  pond  adjacent 
to  the  power  plant,  therefore,  becomes  an  even  larger  relative  contributor  to  ice  fog  than  the  river 
discharge  in  Fairbanks.*  The  cooling  pond  is,  however,  probably  the  least  offensive  of  all  ice  fog 
sources  since  it  contributes  only  ice  crystals.  Products  of  combustion,  which  constitute  the  offen- 
sive portion  of  most  other  sources  of  ice  fog,  are  not  a part  of  the  fog  rising  from  the  pond. 

On  a cold  morning  (-30°C  or  colder)  a cloud  of  ice  fog  rising  from  the  pond  can  be  seen.  It  tends 
to  rise  vertically  since  there  is  rarely  a discernible  wind  during  very  cold  weather.  The  only  air  move- 
ment on  very  cold  days  is  "air  drainage"  flow  (OhUke  1970)  which  generally  moves  from  the  hills 
in  the  east,  across  the  base  toward  the  west.  Since  the  cooling  pond  is  located  on  the  east  edge  of 
the  base,  it  is  in  position  to  contribute  to  the  general  ice  fog  condition.  Furthermore,  it  is  located 
so  that  its  fog  moves  into  the  runway  area. 

The  mechanism  of  formation  of  ice  fog  has  been  studied  at  length  and  is  well  documented. 

Ohtake  gives  the  following  description  of  the  formation  of  ice  fog  (Ohtake  1970,  p.  1 32-1 33): 


♦It  must  be  emphasized  that  producing  heating  steam  does  not  load  the  power  plant;  instead  it 
makes  practical  use  of  energy  that  cannot  be  used  for  electrical  generation,  and  which  would  other- 
wise be  wasted  in  the  cooling  pond,  while  forming  more  ice  fog.  If  enough  “waste"  energy  could 
be  used  for  producing  heating  steam,  the  cooling  pond  would  be  unnecessary. 
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Aerial  photographs  showed  steam  or  water  clouds,  which  are  important  sources  of  ice- 
fog  moisture,  coming  from  the  open  water  of  the  river,  slough  and  cooling  ponds, 
power  plants  and  private  heating  systems.  Such  water  droplets  will  freeze  in  several 
seconds  within  a distance  of  3 to  5 m from  open  water  under  low  temperature  condi- 
tions. ...The  humidity  in  ice  fog  lies  between  water-  and  ice-saturation,  allowing  the 
ice-fog  crystals  to  grow,  but  very  slowly  because  of  the  small  differences  of  saturation 
vapor  pressure  between  them.  This  results  in  ice-fog  crystals  having  the  smallest  size 
of  ice  crystals  and  being  suspended  in  the  air  for  a long  time.  The  modal  size  of  the 
ice  crystals  changes  with  the  changes  in  the  moisture  content  of  the  environmental  air. 

In  an  area  which  has  less  moisture  than  ice  saturation,  the  crystal  sizes  will  be  smaller 
and  smaller  until  the  crystals  disappear.  ...The  most  important  factor  in  the  formation 
of  steam  or  water  droplets  is  not  the  concentration  of  condensation  nuclei  in  the  case 
of  formation  of  steam  from  open  water,  heating  plants  and  car  exhausts  under  ice-fog 
conditions  but  rather  the  temperature  differences  between  water  (not  ice)  and  ambient 
air  temperature. 

A typical  example  of  ice  fog  formation  at  Eielson  AFB  occurred  on  1 2 January  1 973,  when  the 
air  temperature  was  -34°C.  Light  ice  fog  was  evident  along  the  streets  of  the  base,  but  visibility 
there  was  generally  several  miles.  The  air  was  drifting  from  the  north-northeast.  The  fog  cloud 
around  the  pond  was  visible  for  several  miles  as  it  rose  from  the  pond  and  drifted  in  an  ever-widen- 
ing band  across  the  southern  portions  of  the  base.  The  band  was  estimated  to  average  800  m in 
width  and  was  approximately  20  m deep,  as  observed  from  outside  the  cloud.  Upon  entering  the 
cloud  of  fog,  visibility  dropped  to  less  than  50  m,  the  sky  was  obscured,  and  the  sun  was  not  visible. 
(These  limitations  in  visibility  cause  severe  restrictions  in  aircraft  operations  and  make  automobile 
driving  hazardous.)  During  this  period,  the  plume  from  the  power  plant  stacks  was  observed  to 
rise,  then  bend  over,  as  is  typical  during  strong  temperature  inversions  in  the  air  layers  near  the  sur- 
face. The  plume  eventually  rejoined  the  cloud  of  fog  on  the  ground.  This  junction,  however,  was 
well  beyond  the  limits  of  the  base  and  did  not  add  to  fog  on  base. 

An  efficient  approach  to  the  problem  of  ice  fog  suppression  lies  in  suppressing  ice  fog  at  the 
source.  Several  possible  approaches  have  been  proposed.  They  include  the  closed-circuit  cooling 
of  the  condenser  water  and  the  covering  the  cooling  pond  with  different  covers,  including  poly- 
ethylene, styrofoam  balls,  floating  metal  plates,  and  inflatable  buildings.  Injection  wells  have  also 
been  recommended  as  a means  of  disposing  of  cooling  water,  with  new  cooling  water  being  drawn 
from  the  water  table.  Several  of  thv  ' ideas  have  merit  and  should  be  explored  in  depth. 

SECTION  II.  EVAPORATION 

Evaporation  accounts  for  a major  portion  (=*  25%)  of  the  heat  dissipated  from  an  open  water 
surface  during  winter  months.  It  is  also  recognized  as  a significant  contributor  to  the  ice  fog  in  the 
Fairbanks  area.  Elimination  of  this  source  of  ice  fog  by  evaporation  suppression  deprives  the  pond 
of  one  of  the  means  of  dissipating  waste  heat,  so  that  for  the  pond  to  remain  effective,  alternate 
means  of  transferring  this  heat  to  the  atmosphere  or  to  outer  space  must  be  provided.  Therefore, 
a knowledge  of  the  order  of  magnitude  of  evaporative  heat  loss  is  essential  to  ice  fog  suppression, 
so  that  adequate  alternate  heat  transfer  modes  may  be  incorporated  in  the  system. 

Evaporation  is  a function  of  many  interacting  variables.  Water  temperature  and  surface  condi- 
tions influence  the  vapor  flux  leaving  the  surface,  while  wind,  air  temperature,  and  relative  humidity 
influence  the  flux  of  vapor  returning  to  the  water.  The  interaction  of  these  variables  has  made 
mathematical  predictions  of  evaporation  difficult. 

Most  observers  have  noted  that  wind  is  an  important  parameter.  In  fact,  some  equations  have 
been  proposed  with  form  such  that  if  wind  goes  to  zero,  evaporation  ceases.  This,  of  course,  is 
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impossible  since,  even  in  stagnant  air,  diffusion  accounts  for  a net  flux  of  vapor  away  from  the 
surface. 

A still  greater  evaporation  parameter  is  induced  air  movement.  During  ice  fog  conditions  a unique 
situation  develops:  the  air  near  the  ground  (below  50  m)  is  nearly  quiescent,  so  that  conventional 
wind  measurements  show  aero  velocity  (that  is,  air  movements  are  below  the  measurement  thresh- 
old of  conventional  instruments).  The  air  over  the  warm  open  water  surface  of  a cooling  pond  is 
far  from  still,  however.  The  temperature  differences  between  air  and  water  can  be  greater  than 
50°C.  It  is  visually  evident  that  the  air  near  the  water  surface  is  warmed  and  rises  convectively, 
since  vertical  air  currents  can  be  observed  over  the  pond.*  Therefore,  although  meteorological 
measurements  in  the  vicinity  of  a cooling  pond  may  show  zero  wind,  convective  currents  are  usu- 
ally present  and  they  greatly  influence  evaporation. 

Although  extremely  cold  air  can  hold  very  little  moisture,  the  continued  flow  of  new  air  into  the 
pond  area  to  replace  the  air  that  rises  vertically  in  the  convection  plume  provides  a mass  transfer 
mechanism  that  can  carry  off  large  amounts  of  vapor.  In  addition  to  this,  as  Yen  and  Landvatter 
(1970)  have  found,  considerable  supersaturation  of  vapor  in  cold  air  can  exist  above  a warm  water 
source.  It  is  not  surprising,  therefore,  to  find  high  winter  evaporation  rates. 

It  might  also  be  expected  that  equations  developed  for  warmer,  windier  climates  may  underpre- 
dict the  evaporation  during  times  when  meteorological  measurements  indicate  wind  to  be  absent. 


Bowen’s  equation 

To  begin  a discussion  of  the  various  attempts  to  describe  evaporation  in  mathematical  terms,  it 
is  necessary  to  become  acquainted  with  Bowen's  ratio,  which  is  used  in  the  derivation  of  virtually 
all  modern  evaporation  equations. 

Bowen  (1926)  concluded  that  "the  process  of  evaporation  and  diffusion  of  water  vapor  from  any 
water  surface  into  the  body  of  air  above  it  is  exactly  similar  to  that  of  the  conduction  or  ‘diffusion’ 
of  specific  heat  energy  from  the  water  surface  into  the  same  body  of  air.”  He  proposed  a ratio  RB, 
which  has  become  known  as  Bowen’s  ratio. 

Rb  is  defined  as  the  ratio  of  convective  heat  loss  (?H  to  evaporative  heat  loss  QE.  It  is  assumed 
in  Bowen’s  derivation  that  the  eddy  diffusivities  of  water  vapor  Kw  and  heat  /Ch  are  equal  in  air 
(a  result  derivable  from  kinetic  theory).  Several  investigators  have  reported  satisfactory  results 
using  Bowen’s  ratio  under  stable  and  neutral  conditions  (Pasquill  1948,  Dingman  et  al.  1967,  Ander- 
son 1952).  Pasquill  (1948),  on  the  other  hand,  has  reported  that  under  unstable  conditions,  i.e. 
when  the  water  temperature  is  greater  than  the  overlying  air  temperature  (Tw  ,>  7"a),  Kh  can  be  as 
much  as  three  times  as  large  as  K^„  However,  Brutsaert  (1965)  cites  two  reports  that  concern 
themselves  with  the  subject:  one  which  found  /Ch  = /Cw  regardless  of  stability  conditions,  and  the 
other  which  confirmed  that  they  were  equal  to  each  other  except  under  highly  stable  conditions. 

Working  under  winter  conditions,  Rimsha  and  Donchenko  (1958)  have  successfully  used  Bowen’s 
equation  (with  Kh  = /Cw)  to  calculate  heat  losses  from  ice-free  stretches  of  rivers.  In  view  of  these 
studies  it  seems  reasonable,  for  the  conditions  of  this  study,  to  assume  (as  did  Rimsha  and  Don- 
chenko) that  Bowen’s  ratio  can  be  used  with  Kh  = Kw.  Bowen’s  ratio  ( RB ) is 
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♦The  situation  may  be  analogous  to  the  generation  of  thunderstorms  in  the  summer,  where  mois- 
ture-laden air,  heated  by  a warm  surface,  rises  to  elevations  of  1 3,000  m (40,000  ft)  and  higher. 
Currents  and  turbulence  within  these  columns  of  rising  air  have  been  known  to  destroy  aircraft 
caught  within  them. 


Temperature  ( C) 

Figure  2.  Saturation  water  vapor  pressure  vs  temperature 
(after  Eagieson  1972). 

where  P = atmospheric  pressure  (N  m'2 ) 
rw  = temperature  of  the  water  (°C) 

7,  = temperature  of  the  air  (°C) 

ew  = saturated  vapor  pressure  at  the  temperature  of  the  water  (N  m'2 ) 
et  = vapor  pressure  at  the  temperature  of  the  air  (N  m*2)  (Fig.  2) 

<?H  = convective  heat  loss  (W  m‘2 ) 

<?E  = evaporative  heat  loss  (W  m"2 ) 

CB  = Bowen's  constant  (°C‘I ). 

Bowen  (1926)  evaluated  CB  for  three  wind  conditions.  Case  I represents  a condition  where,  in 
Bowen’s  words,  "the  whole  quantity  of  air  under  consideration  is  completely  changed  in  tempera- 
ture and  moisture  content  to  that  of  the  layer  of  air  in  contact  with  the  water.”  In  this  case,  diffu- 
sion is  not  the  limiting  factor  and,  therefore,  Kh  and  ACW  do  not  enter  into  the  equation.  For  Case  I, 
CB  is  found  empirically  to  be 

CB  = 0.501  °C-‘ . 

Case  II  represents  the  situation  where  diffusion  dominates  completely.  Heat  and  water  vapor  are 
immediately  carried  away  after  diffusing  through  a stationary  layer  in  contact  with  the  water.  For 
this  case  CB  = f(ACh//Cw)  and  was  empirically  determined  to  be 

CB  = 0 442  °C*' . 

Case  III  represents  conditions  where  wind  u is  a function  of  the  height  above  the  water  u = f{y)  ~ 

Ky ",  where  AC  is  a constant.  This  is  the  condition  that  Bowen  considers  the  “average  wind  condi- 
tion” and  for  this  case,  his  analysis  shows  that 

CB  =0.46°C*‘. 


Bowen  concluded  that  the  best  value  for  /?B  utilized  CB  = 0.46.  Using  this  value,  Bowen's  ratio 
Rb  becomes 

, (rw-r) 

= (®-1  x 10”* ) (P)  7~  • (2) 
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Bowen  states  that  this  equation  is  valid  f<#thc  situation  where  Tw  is  low  enough  so  that  the  vol- 
ume of  air  is  not  changed  significantly  by  the  water  vapor  evaporated  into  it,  and  where  condensed 
evaporated  water  does  not  precipitate  back  into  the  pond.  These  two  restrictions  are  both  violated 
to  some  degree  at  a cooling  pond  during  ice  fog  conditions.  This  may  reflect  on  the  accuracy  with 
which  equations  based  on  Bowen's  ratio  can  describe  the  actual  evaporation  at  these  times. 

Evaporation  equation 

Traditionally,  there  have  been  two  basic  approaches  to  develop  the  theory  for  predicting  evapora- 
tion: the  “aerodynamic"  method  and  the  “energy-budget”  method. 

The  aerodynamic  method  approach.  Most  empirically  derived  equations  of  the  aerodynamic 
method  have  the  same  basic  form  as  Dalton's  Law: 

E=\f(u)\(e  w-eJ  (3) 

where  E = the  evaporation  in  a unit  of  time  and  area 

eaj  = the  water  vapor  pressure  in  the  air  at  the  2-m  level 
ew  = the  saturation  vapor  pressure  at  the  temperature  of  water  surface 
f(u)  = a function  of  the  windspeed. 

The  form  of  the  wind  function  is  usually  taken  such  that 
f(u)  = (a+bu) 

where  a and  b arc  empirically  derived  constants  and  u is  the  average  windspeed  2 m above  the 
surface. 

At  ice  fog  temperatures  e ^ can  be  replaced  by  the  saturated  vapor  pressure  ea  without  significant 
error,*  and  therefore  eq  3 can  be  rewritten: 


E = (a+bu)  (ew-ea)  (4) 

where  ea  = saturated  vapor  pressure  at  the  temperature  of  the  air.  The  constants  a and  b must  be 
evaluated  experimentally. 

Dingman  et  al.  (1967)  propose  a form  of  Dalton’s  equation  which  they  attribute  to  Kohler’s 
work  that  is  “a  simple,  well-accepted  formula  for  calculating  daily  evaporation  rates”: 

E = (0.525x  1 0‘J  + 1 ,0007x  1 0'J u)  (ew-ea ) (5) 

where  u = windspeed  in  m s'1  and  ew  and  ea  are  in  millibars. 


♦It  must  be  remembered  that  this  approximation  will  only  be  usable  at  temperatures  well  below 
0°C.  For  water  temperatures  of  10°C,  the  maximum  error  introduced  by  this  approximation  is 
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Vapor  pressure 
at  water  terrip. 
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Vapor  pressure 
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-20 

12.3 

1.25 

11.05 

12.3 

10 

-30 

12.3 

0.51 

11.79 

12.3 

4 

-40 

12.3 

0.19 

12.11 

12.3 

1.5 

Maximum  errors  will  be  produced  if  the  air  entering  the  pond  area  is  completely  dry  (eaj  = 0),  and  j 

air  leaving  the  pond  is  saturated  (relative  humidity  = 1 00%).  If  the  air  entering  is  close  to  satura-  ( 

tion,  the  error  approaches  zero. 
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Multiplying  by  the  latent  heat  of  vaporization,  eq  S becomes  an  expression  of  the  evaporation 
heat  loss  Qr : 

<?E  = (1.51 8 + 2.91 1t/)(ew-ea)  (6) 

where  QE  = evaporation  heat  loss  (W  m'J ) 
u = wind  velocity  (m  s'1 ) 

ew  = saturation  vapor  pressure  at  water  temperature  (mb) 
e,  = saturation  vapor  pressure  at  air  temperature  (mb). 

Using  Dalton’s  equation  and  evaluating  the  constants'!?  and  b,  Behlke  and  McDougall  (1973) 
found  that  values  of  a = 1 3.82  and  b = 0 gave  the  optimum  correlation  to  their  measured  data. 

Their  results  in  the  form  of  evaporation  heat  loss  rate  were 

(?E  = 13.82(ew-ea)  (7) 

where  (?E  is  in  W m'J,  ew  is  in  mb  and  ea  is  in  mb. 

Their  study  is  of  particular  interest  since  it  was  conducted  in  interior  Alaska  during  one  of  the 
same  winters  that  the  cooling  pond  studies  (of  this  paper)  at  Eielson  AFB  were  in  progress.  The 
two  investigations,  however,  differed  markedly.  The  Behlke-McDougall  experiments  were  con- 
ducted on  the  roof  of  a four-story  building  located  on  a hill  above  the  Fairbanks  valley  floor,  and 
were  influenced  by  the  surrounding  buildings  and  the  local  "heat  island”  of  the  University  of  Alaska 
complex.  The  evaporation  studies  at  Eielson  were,  in  contrast,  isolated  from  all  other  structi,rcs  by 
at  feast  300  m (1000  ft).  This  site  was  on  the  valley  floor  and,  additionally,  the  evaporation  pans 
used  in  this  experiment  were  subject  to  the  microclimate  of  the  open  water  of  the  pond.  Neverthe- 
less, it  is  interesting  to  note  that  in  spite  of  these  differences,  neither  study  showed  appreciable 
wind  effect  (p.  16).  The  conclusion  follows  that,  during  ice  fog  conditions  in  the  Fairbanks  region, 
the  wind  is  not  a significant  factor  in  evaporation  since  it  is  absent  during  these  times.  If  wind  were 
present,  it  would  not  only  cause  even  greater  evaporation  but  would  also  disperse  the  ice  fog,  pos- 
sibly alleviating  this  problem. 

Another  study  using  the  Dalton  equation  approach  was  reported  by  two  Russian  scientists  in- 
vestigating ice  formation.  Rimsha  and  Donchenko  (1958)  concluded  (after  comparisons  of  actual 
heat  loss)  that  for  winter  conditions  on  rivers  and  reservoirs  throughout  Russia,  the  convective  heat 
loss  rate  (?H  could  be  expressed  by 


<?H=(^  + >-89")(7w-7a)  (8) 

where  u = the  windspeed  in  m s'1  at  2 m above  the  surface  and  Fc  = the  portion  of  heat  trans- 
ferred by  free  convection.  Fc  was  determined  empirically  by  them  to  be 

Fc  = 3.87  +0.1 7(TW-Ta).  (9) 

Applying  the  Bowen  ratio  (assuming  pressure  = 1000  mb),  the  rate  of  evaporative  heat  loss  Qv  is 
determined  to  be 


Qr : = (0J5SFc+2.94u)(ew-eJ.  (10) 

Combining  eq  9 and  10  gives 

<?F  ={6.04+0.264(7w-ra)  + 2.944/|(ew-ea).  (11) 
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This  is  often  called  Rimsha  and  Donchenko’s  “Russian  Winter  Equation.”  The  advantage  of  this 
equation  to  the  needs  of  this  investigation  lies  in  the  term  FQ.  During  ice  fog  conditions  when  wind 
is  absent,  eq  10  reduces  to 

\ - (0.755Fc)(ew-ea).  (12) 

Evaporation  is  thus  expressed  as  a function  of  convective  conditions  represented  by  Fc.  Since 
F(  was  determined  for  winter  conditions  in  Russia,  the  applicability  of  Fc  in  the  revised  Rimsha 
and  Donchenko  equation  (eq  1 2)  may  represent  conditions  at  the  Eielson  Pond  more  closely  (in 
concept  at  least)  than  any  of  the  other  equations  discussed.  Their  studies  were  addressed  to  open 
water  on  rivers  and  reservoirs;  however  (and  this  may  be  a significant  difference),  ice-free  areas  on 
rivers  and  reservoirs  expose  water  that  may  be  warmer  *han  the  surrounding  air,  but  many  degrees 
colder  than  the  water  surface  of  a cooling  pond  (12.6°C  was  the  average  winter  temperature  during 
these  experiments  for  the  Eielson  pond).  In  addition,  the  Russian  studies  were  conducted  at  sites 
throughout  Russia,  many  of  which  (e.g.  Chu  River,  Angara  River)  are  too  far  south  to  be  considered 
arctic.  Thus  the  Russian  Winter  Equation  might  be  expected  to  yield  somewhat  different  values  than 
the  Eielson  studies. 

The  energy  budget  approach.  The  energy  budget  approach  to  evaporation  measurements  is 
merely  an  application  of  the  concept  of  conservation  of  energy.  Williams  (1963,  p.  32)  quotes 
Penman,  “The  fundamental  basis  of  the  energy  budget  approach  is  unchallenged;  the  challenge  is 
to  our  ability  to  measure  or  estimate  all  the  quantities  needed  to  exploit  the  principles  of  the  con- 
servation of  energy.” 

In  mathematical  terms  the  method,  when  applied  *o  a cooling  pond,  can  be  stated 

^in+^sw  + ^o«lier-^l»”^E-<?H  = 03) 

where  Qjn  = rale  of  cooling  water  heat  energy  input  to  the  pond 

(?sw  = net  shortwave  radiation  rate  incoming  from  the  sun  (direct  and  diffuse) 

Mother  = a*'  olhcr  energy  input  rates  not  considered  in  the  first  two  terms  (i.e.  seepage,  ground 
water  make-up,  precipitation,  outflows,  etc.) 

<?.w  = net  longwave  radiant  heat  loss  rate 

Qt  = evaporation  heat  loss  rate  (convccted  and  conducted  from  the  systems) 

Qh  = convective  heat  loss  rate 

<?,  = rate  of  increase  of  heat  energy  stored  in  the  body. 

It  is  immediately  obvious  that  this  method  is  dependent  on  the  accuracy  of  a number  of  com- 
ponents, some  of  which  arc  difficult  if  not  impossible  to  evaluate  (e.g.  seepage  of  warm  springs 
into  the  bottom  of  the  pond).  In  addition,  because  of  the  problems  in  accurately  assessing  the 
change  in  the  heat  storage  term  Qs,  this  method  is  difficult  to  apply  accurately  (especially  for  short 
periods  of  time).  This  is  particularly  true  for  small  bodies  of  water  such  as  cooling  ponds,  since  the 
heat  storage  can  and  does  change  rapidly,  due  to  fluctuations  of  both  pond  level  id  heat  input 
from  the  incoming  cooling  water.  However,  this  method  has  been  used  with  apparent  success  in 
studies  at  Lake  Hefner  for  time  periods  of  seven  days  or  greater  (Kohler  et  al.  1955).  The  method 
would  probably  also  be  adequate  when  applied  to  cooling  ponds  if  time  periods  were  long  enough 
to  average  short-term  storage  changes. 

In  an  effort  to  derive  a relationship  that  is  based  upon  a limited  number  of  reasonably  easy-to- 
obtain  meteorological  values,  Penman  (1948)  combined  the  energy  balance  approach  with  the  aero- 
dynamic approach.  The  derivation  is  lengthy  and  the  reader  is  referred  to  Penman’s  paper  listed  in 
the  bibliography.  The  Penman  equation  is  usually  written  as 


(14) 


°*‘I?7(0'4n£-1) 


where  Qe  = evaporative  heat  loss  rate  (W  m*2 ) 

A = slope  of  the  curve  of  saturation  vapor  pressure  (mb)  vs  temperature  (°C);  [d(e)/dT]  at 
air  temperature  Ts 

Et  = the  evaporation  given  by  an  equation  of  the  form  of  eq  4 (W  m*2 ) 

Qt  = net  radiant  energy,  expressed  in  the  units  of  EtL 

1 »«b[M,)/M.)I  =cBp  (mb>  (°c'') 

L * latent  heat  of  evaporation  (cal  g'1 ). 

An  advantage  Penman  claims  for  this  approach  is  that  “it  eliminates  the  need  for  water  surface 
temperature,  the  parameter  most  difficult  to  measure”  (Penman  1 948,  p.  1 20).  However,  the 
method  is  based  on  the  assumption  that  heat  storage  of  the  system  does  not  change.  Therefore,  it 
must  be  applied  over  longer  periods  of  time  in  order  to  “average  out”  the  effects  of  storage  fluctua- 
tions. 

The  term  y is  a function  of  the  eddy  diffusivity  ratio  Kh/Kw  discussed  earlier,  and  if  this  ratio 
departs  from  the  assumed  value  of  1 , adjustments  of  y must  be  made  to  compensate.  Another  dis- 
advantage of  the  method  is  that  it  relies  on  the  standard  Dalton  relationship  for  the  wind  function, 
and  values  for<7  and  b in  eq  4 must  be  evaluated.  It  does  not  consider,  in  its  present  state,  the  spe- 
cial case  of  convective  instability. 


Evaporation  measurements  during  ice  fog  conditions 

Very  few  experiments  have  been  conducted  which  measured  evaporation  under  actual  ice  fog 
conditions.  Ohtake  (1970)  measured  the  evaporation  of  water  from  two  shallow  pans  (13.4  and 
30.5  cm  in  diameter  and  approximately  1 cm  deep).  The  pans  were  exposed  to  the  atmosphere 
and  kept  from  freezing  with  an  infrared  lamp.  Water  temperature  in  the  pans  was  measured  with 
thermistors,  and  was  maintained  above  freezing  by  adjusting  the  distance  from  the  lamp  to  the  pan. 
Measurements  were  made  at  air  temperatures  ranging  from  -18°  to  -40°C.  Water  temperature  was 
“...kept  around  4°  to  10°C  with  an  occasional  maximum  of  15°C...”  (Ohtake  1970,  p.  118).  A 
total  of  10  measurements  yielded  an  average  evaporation  of  5000  g m*2  day’1 . The  rather  large 
variations  of  water  temperatures,  coupled  with  the  short  time  increments,  gave  rather  wide  scatter 
in  evaporation  measurements,  which  ranged  from  2380  to  8530  g m*2  day*1 . In  contrast  to  the 
small  scale  of  Ohtake’s  experiment,  a cooling  pond  is  large  enough  to  create  its  own  microclimate, 
and  enough  heat  may  be  dissipated  from  the  large  open  water  surface  of  the  pond  to  create  sub- 
stantial convective  currents.  Such  a microclimate  is  not  nearly  as  well  developed  in  the  case  of 
small  pans  on  land,  and  the  difference  suggests  that  the  values  may  not  be  completely  representa- 
tive of  conditions  in  the  larger  system.  Anderson  (1952)  and  Kohler  ct  al.  (1955)  discuss  this 
point  at  length. 

Yen  and  Landvatter  (1970)  measured  the  evaporation  from  stretched  blotter  paper  with  below- 
freezing  air  flowing  over  it.  Their  measurements  covered  a water  temperature  range  from  0.8°C  to 
30.4°C,  and  air  velocities  from  0.10  to  1 .60  m s*1 . Their  apparatus  allowed  the  variables  involved 
in  evaporation  to  be  controlled  in  a much  better  fashion  than  is  available  in  field  experiments.  A 
total  of  64  runs  were  reported.  Evaporation  ranged  from  1296  to  7058  g m"2  day-1  for  water 
temperatures  of  1 .7  f to  6.7°C  and  below-freezing  air,  and  from  5642  to  24,019  g m*2  day*’  for 
water  temperatures  between  24.9°  and  30.3°C  and  above-freezing  air  temperatures.  Unfortunately, 
no  measurements  were  made  for  water  temperatures  between  8.3°  and  24.9°C,  the  range  in  which 
the  Eielson  cooling  pond  water  surface  temperatures  fall  most  of  the  time.  The  emphasis  in  these 
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Figure  3.  Colorado  evaporation  pans. 


experiments  on  forced  convection  suggests  that  they  may  not  be  closely  representative  of  the  con- 
ditions on  the  cooling  pond  where  air  movement  is  predominantly  free  convection.  The  scale  in 
this  case,  also,  is  different.  Nevertheless,  these  measurements  represent  evaporation  under  condi- 
tions which  are  more  nearly  similar  to  those  predominating  at  the  cooling  pond  than  for  most  of 
the  evaporation  measurements  reported  in  the  literature,  and  agree,  in  order  of  magnitude  at  least, 
to  those  which  were  taken  at  the  Eielson  cooling  pond.  Table  I points  this  out.  If  only  those  runs 
with  below-freezing  air  are  considered,  the  average  values  of  Ohtake,  and  Yen  and  Landvatter  are 
within  28%  of  each  other,  with  Yen  and  Landvatter  showing  the  lower  average.  On  the  other  hand, 
if  the  higher  water  temperature  and  above-freezing  air  values  of  Yen  and  Landvatter’s  study  are  in- 
cluded, then  their  average  value  is  greater  than  Ohtake’s  by  27%.  This  suggests  that  good  agree- 
ment of  the  values  of  evaporation  might  have  been  obtained  if  both  studies  had  covered  the  same 
temperature  range. 

Evaporation  measurements  at  Eielson  Power  Plant  cooling  pond 

Two  standard  Colorado  floating  pans  were  placed  in  the  pond  (Fig.  3).  The  pans  were  read  con- 
currently in  the  same  location  and  the  values  were  compared  (Table  I).  The  average  difference  be- 
tween all  values  for  the  two  pans  increased  as  spring  approached.  In  January  an  average  difference 
of  2.1%  was  recorded  for  measurements  covering  a period  of  1 1 8 hours,  but  in  March  the  difference 
between  the  pans  had  increased  to  1 7.6%  over  a period  of  1 66  hours.  This  is  probably  a function  of 
a number  of  variables.  Air  movement  across  the  pond,  increased  solar  radiation  as  spring  approached, 
and  warm  currents  in  the  pond  could  all  have  contributed  to  the  difference.  As  solar  radiation  in- 
creased, the  degree  to  which  the  pans  were  in  the  sun  would  influence  their  evaporation.  If  the 
shadow  from  the  pump  house  on  the  pond  (Fig.  1)  selectively  shaded  one  pan  and  not  the  other,  a 
differential  evaporation  would  result.  Another  possibility  is  that  warm  currents  developed  in  the 
pond  (from  the  nearby  warm  water  diversion  inlet);  should  one  pan  be  selectively  warmed  by  such 
a current,  it  would  experience  higher  evaporation.  Finally,  when  the  pans  were  in  close  proximity 
to  each  other,  air  movement  from  either  the  east  or  west  could  affect  the  first  pan  it  encountered; 
at  the  same  time,  this  pan  would  provide  a shelter  for  the  second  pan.  As  weather  warmed  from 
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Table  I.  Evaporation  data  from  standard  Colorado  floating  pans, 
and  smaller  pans  prepared  for  this  study. 


Date 

Pan 

no. 

Duration 

(hr) 

Air 

temp. 

CO  _ 

Water 

temp. 

Cc) 

A vg  wind 
velocity 
(ms'1) 

Evaporation 
(g  nr2  day'  ) 

Colorado  floating  pan  data 

2 Jan  73 

2 

95.0 

-16.0 

8.4 

0 

3720 

5 Jan  73 

2 

72.5 

-16.0 

11.1 

0 

3812 

11  Jan  73 

1 

23.5 

-11.0* 

10.8 

0 

3633 

2 

25.0 

-11.0* 

10.9 

0 

4079 

12  Jan  73 

1 

25.1 

-31.0 

9.5 

0 

3996 

2 

24.8 

-31.0 

9.0 

0 

3972 

13  Jar.  73 

2 

24.2 

-38.5 

10.2 

0 

6853 

14  Jan  73 

1 

25.2 

-38.0 

9.0 

0 

5485 

2 

24.2 

-38.0 

9.3 

0 

5622 

15  Jan  73 

1 

23.5 

-40.0 

8.4 

0 

5728 

2 

23.8 

-40.0 

8.6 

0 

5476 

18  Jan  73 

1 

69.8 

-33.0 

9.5 

0 

5247 

2 

70.0 

-33.0 

9.3 

0 

5013 

20  Jan  73 

1 

46.8 

-36.0 

7.5 

0 

6738 

2 

48.2 

-36.0 

8.8 

0 

6603 

21  Jan  73 

1 

28.5 

-39.0 

8.2 

0 

4854 

23  Jan  73 

1 

46.8 

-38.0 

8.2 

0-0.50 

7656 

2 

49.5 

-38.0 

8.5 

0-0.50 

3980 

6 Feb  73 

1 

25.0 

-18.0 

10.0 

0.20 

3364 

2 

24.5 

-18.0 

10.0 

0.20 

3731 

7 Feb  73 

1 

24.0 

-18.0 

12.2 

0.05 

4676 

2 

24.5 

-18.0 

12.0 

0.05 

5617 

8 Feb  73 

1 

23.8 

-17.0 

13.4 

0.15 

7700 

2 

24.0 

-17.0 

13.4 

0.15 

4676 

9 Feb  73 

1 

23.6 

-20.0 

13.5 

0 

2821 

2 

?3.5 

-20.0 

13.5 

0 

7488 

10  Feb  73 

1 

26.7 

-15.0 

12.5 

0 

8301 

2 

26.2 

-15.0 

12.5 

0 

5894 

11  Feb  73 

1 

25.2 

-19.0 

11.0 

0.90 

4440 

2 

25.2 

-19.0 

11.0 

0.90 

7508 

12  Feb  73 

1 

24.0 

-21.0 

10.5 

3.00 

5495 

13  Feb  73 

2 

22.8 

-26.0 

11.5 

1.50 

5823 

15  Feb  73 

1 

49.5 

-22.0 

13.5 

0.12 

6253 

2 

49.5 

-22.0 

' 14.0 

0.12 

£184 

21  Feb  73 

1 

24.0 

- 5.0 

15.6 

1.40 

6105 

2 

23.8 

- 5.0 

15.7 

1.40 

4834 

22  Feb  73 

1 

24.8 

- 4.0 

14.9 

0.10 

5393 

2 

24.8 

- 4.0 

15.0 

0.10 

5996 

23  Feb  73 

1 

23.5 

- 4.5 

13.8 

0 

4183 

2 

23.8 

- 4.5 

14.0 

0 

4990 

24  Feb  73 

1 

23.0 

-12.0 

13.6 

0.62 

5141 

2 

22.8 

-12.0 

13.6 

0.62 

3835 

26  Feb  73 

1 

48.7 

-24.0 

11.0 

0.50 

4669 

2 

49.0 

-24.0 

11.1 

0.70 

4604 

♦These  temperatures  are  from  Air  Force  measurements  at  Eielson  AFB,  approximately  2 miles 
from  the  pond  site. 
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Table  I (corn’d). 


Date 

Pan 

no. 

Duration 

(hr) 

Air 
temp . 

rc) 

Water 

temp. 

rc) 

Avg  wind 

velocity 

(ms-') 

Evaporation 
( gm day-'1) 

27  Feb  73 

1 

22.8 

-16.0 

5.4 

0.70 

5835 

2 

22.4 

-16.0 

5.3 

0.70 

5971 

3 Mar  73 

2 

23.2 

-11.0 

7.0 

0.35 

1931 

8 Mar  73 

1 

24.0 

-10.0 

9.7 

0.20 

3839 

2 

24.5 

-10.0 

10.0 

0.20 

3732 

9 Mar  73 

1 

26.0 

- 9.5 

12.4 

0.30 

1965 

2 

25.6 

- 9.5 

12.6 

0.30 

2603 

10  Mar  73 

1 

22.3 

-13.0 

12.6 

1.60 

4270 

2 

22.8 

-13.0 

12.9 

1.00 

4725 

1 1 Mar  73 

1 

24.20 

- 9.5 

13.3 

0.15 

5420 

2 

23.75 

- 9.5 

13.4 

0.15 

4363 

12  Mar  73 

1 

22.20 

-15.0 

13.6 

1.15 

5740 

2 

22.60 

-15.0 

13.9 

1.15 

6672 

14  Mar  73 

1 

50.00 

-12.0 

14.3 

0.70 

6045 

2 

49.50 

-12.0 

14.4 

0.70 

5856 

Small  pan  data 

27  Jan  72 

S-1 

6 

18 

10.2 

1.50 

2850 

S-2 

6 

18 

10.2 

1.50 

3270 

28  Jan  72 

S-1 

4 

16 

10.2 

0 

2440 

S-2 

4 

16 

10.2 

0 

2820 

1 Feb  72 

S-1 

5 

14 

0.5 

1.50 

1122 

S-2 

5 

14 

0.5 

1.50 

1230 

1 Mar  72 

S-1 

4 

-23 

3.0 

0 

2119 

S-2 

4 

-23 

3.0 

0 

2650 

S-3 

4 

-23 

3.0 

0 

2119 

2 Mar  72 

S-1 

6 

-26 

2.0 

0 

3470 

S-2 

6 

-26 

2.0 

0 

3060 

3 Mar  72 

S-1 

3 

-24 

2.0 

0 

1410 

S-2 

3 

-24 

2.0 

0 

3760 

S-3 

3 

-24 

2.0 

0 

3290 

S-4 

3 

-24 

2.0 

0 

3530 

7 Mar  72 

S-1 

4 

-20 

5.0 

0.75 

7900 

S-2 

4 

-20 

4.0 

0.75 

4690 

S-3 

4 

-20 

4.0 

0.75 

3670 

S-4 

4 

-20 

2.0 

0.75 

5770 

8 Mar  72 

S-1 

4 

-14 

4.5 

3.00 

3840 

S-2 

4 

-14 

3.0 

3.00 

3770 

10  Mar  72 

S-1 

5 

-20 

5.5 

0 

4000 

S-2 

5 

-20 

5.0 

0 

3300 

S-3 

5 

-20 

4.5 

0 

1890 

14  Mar  72 

S-1 

5 

-18 

6.5 

0 

4070 

S-2 

5 

-18 

6.0 

0 

4240 

S-3 

5 

-18 

5.5 

0 

3110 

S-4 

5 

-18 

5.5 

0 

4240 

\ 
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Table  II.  Evaporation  rates  during  ice  fog  temperatures. 


I 


i 


Average 

evaporation 

(g  m'2  day'1 ) (mm  day 1 ) 

Water  temp, 
range 

rc) 

No. 

of 

meas. 

Data 
spread 
(mm  day') 

Ohtake  (1970) 

5040 

5.0 

4 -15 

10 

0.9-  5.9 

Behlke  and  McDougall  (1973) 

4464 

4.5 

2.7-11.9 

10 

2.6-  4.3 

Yen  and  Landvatter  (1970)' 

6912 

6.9 

0.8-30.3 

64 

0.9-16.7 

Yen  and  Landvatter  (1970)6 

5698 

5.7 

0.8-30.3 

— 

- 

Yen  and  Landvatter  (1970)2 

3600 

3.6 

0.8-  7.9 

37 

0.9-  4.9 

Eielson  cooling  pond  studies3 

5040 

5.0 

0.0-16.0 

79 

1.3-  5.8 

Russian  winter  equation4 

5393 

5.4 

10 

- 

- 

Revised  Russian  equation4 

4314 

4.3 

10 

- 

- 

Eielsor  form,  Dalton  equation4 

5045 

5.0 

10 

— 

- 

Kohler  equation4 

585 

0.5 

10 

- 

- 

Modified  Penman  equation4 

3456 

3.4 

calculated 
for  12° 

— 

Thompson  (1970)5 

2736 

2.7 

calculated 

- 

- 

1 Laboratory-  measurements,  all  values  included. 

2Only  values  for  air  temperatures  below  freezing  included. 

3 Values  for  the  Eielson  cooling  pond  were  the  average  values  for  all  conditions  during  the  testing 
throughout  the  winter. 

4 Values  were  calculated  for  -20 °C  air  temperatuie  and  +10 °C  water  temperature. 

5 Using  Thomwaite  and  Holzman’s  equation  (reported  by  Ohtake  1970). 

6 Composite  average  of  above-freezing  and  below-freezing  air  temperatures  with  water  tempera- 
tures below  8.4°C. 


January  to  March,  both  the  magnitude  and  frequency  of  wind  increased  over  the  pond  (Table  I), 
which  may  help  explain  the  increase  in  the  difference  between  the  pans  as  spring  approached.  The 
overall  difference  between  the  pans  of  less  than  9%,  however,  is  within  reasonable  expectations  for 
evaporation  field  measurements  even  under  temperate  conditions  (see  data  spread  for  laboratory 
measurements  of  Yen  and  Landvatter,  Table  II).  The  total  spread  in  evaporation  measurements  at 
Eielson  was  less  than  reported  by  either  Yen  and  Landvatter  (1970),  or  Ohtake  (1970)  in  their 
measurements  (see  Table  II). 

The  rate  of  evaporation  from  the  pans  was  determined  by  measuring  the  volume  of  water  evapo- 
rated each  day.  This  was  done  by  raising  the  pan  partially  out  of  the  water  with  a lift,  and  leveling 
it.  Water  was  then  added  to  the  level  of  a point  gage  installed  in  the  pan.  The  added  volume  of 
water  was  measured  with  a graduated  cylinder  to  the  nearest  10  ml.  Precision  of  the  point  gage  in 
the  pans  was  determined  experimentally  in  the  laboratory  and  found  to  have  a standard  deviation 
of  38  ml  over  40  measurements.  Evaporation  volumes  for  each  day  were  in  the  order  of  several 
thousand  ml  in  each  pan,  so  that  measurement  errors  were  less  than  1%. 

To  obtain  data  on  evaporation  from  cooler  surface  areas  of  the  pond,  four  33-cm-diam  pans 
were  used  to  measure  evaporation  at  different  locations  on  the  pond  (Fig.  4).  These  smaller,  more 
mobile  pans  were  compared  with  each  other  and  with  the  larger  Colorado  pans.  Measurements  were 
taken  in  the  same  manner  as  the  Colorado  pans,  except  that  the  period  of  measurement  was  much 
shorter  since  the  pans  were  in  colder  water  and  would  freeze  when  left  out  overnight.  It  was  not 
convenient  to  stabilize  the  smaller  pans  before  measurement;  therefore,  the  accuracy  of  these  data 
points  is  subject  to  more  influence  from  outside  factors  (such  as  waves  which  rocked  the  pan 
slightly,  causing  reflections  that  made  measurement  difficult). 
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Figure  4.  Small  evaporation  pans. 


Fifty-seven  measurements  were  made  with  the  Colorado  pans,  covering  a time  span  of  495  hours 
of  evaporation.  The  smaller  pans  were  used  over  short  time  spans,  usually  not  exceeding  6 hours. 
Since  the  temperature  of  the  water  in  the  small  pans  often  approached  0°C,  ice  sometimes  formed 
over  a portion  of  the  surface.  This  was  stirred  away  at  the  hourly  checks.  The  average  water  temp- 
erature in  these  pans  was  3.8°C,  while  the  large  Colorado  pans,  positioned  in  a warmer  part  of  the 
pond,  averaged  1 1°C. 

Traditionally,  a "pan  constant”  for  the  evaporation  pans  has  been  used  to  account  for  the  differ- 
ence in  temperature  between  the  water  in  a lake  or  river  and  that  in  the  pan.  Kohler  (1954,  p.  139) 
states,  “The  practice  of  applying  a coefficient  to  pan  evaporation  to  estimate  that  occurring  from  an 
adjacent  lake  is  of  long  standing,  even  though  the  reliability  of  this  approach  has  always  been  sub- 
ject to  question.”  Since  the  water  temperature  in  the  pans  under  measurement  conditions  encoun- 
tered during  most  studies  was  warmer  than  that  of  the  adjoining  lake,  a pan  coefficient  of  less  than  1 
is  applied  (values  of  0.86  to  0.95  are  usually  used  for  floating  Colorado  pans)  (Kohler  1954,  Ander- 
son 1952).  In  the  measurements  made  at  Eielson  AFB,  however,  the  pan  water  temperature  was 
colder  than  that  of  the  surrounding  pond  in  all  cases  (averaging  2.47°C  lower). 

This  would  indicate  the  use  of  a pan  constant  greater  than  1 . However,  since  the  temperature 
over  the  surface  of  the  pond  graded  from  0°  at  the  edge  of  the  ice  cover  to  as  high  as  16°C  near 
the  inlet,  the  pans  can  be  assumed  to  indicate  evaporation  from  those  areas  of  the  pond  where  the 
surface  temperature  is  the  same  as  that  of  the  pans.  Since  the  average  pan  temperature  is  slightly 
higher  t'  an  the  average  temperature  of  the  open  water  surface  of  the  pond,  a pan  constant  of  less 
than  1 might  be  indicated. 

Another  factor  should  be  considered.  The  surface  of  the  water  in  the  pond  lies  approximately 
2 m below  the  dikes  that  make  up  the  sides.  As  air  moves  across  the  pond,  the  dikes  produce  a 
sheltered  area  on  the  lee  side  (that  is  approximately  10  m wide).  This  area  is  affected  very  little 
by  the  wind,  and  is  visible  as  still  water  next  to  rough  water,  even  during  light  breezes.  Although 
there  weren’t  even  light  breezes  during  ice  fog  events,  a sizable  portion  of  the  evaporation 


data  was  taken  during  warmer  weather  wnen  breezes  were  present.  The  remaining  part  of  the  pond 
(70  m)  received  the  effect  of  the  wind.  The  floating  evaporation  pans  have  a similar  configuration 
as  the  pond  with  respect  to  the  wind  at  their  sides.  The  sides  of  the  pans  extend  approximately 
10  cm  above  the  pond  surface  and  shelter  the  surface  of  the  enclosed  water,  in  much  the  same  man- 
ner as  the  dikes  shelter  the  pond.  The  pans,  however,  have  a somewhat  different  geometry  than  the 
pond,  and  therefore,  the  effect  of  the  wind  on  evaporation  from  the  pans  is  most  likely  a bit  less 
than  that  on  the  pond  itself.  This  would  indicate  the  use  of  a pan  constant  slightly  greater  than  1, 
at  least  for  windy  periods. 

The  horizontal  air  movement  during  ice  fog  events  is  so  slight,  however,  that  the  shadow  effect 
of  the  dike  was  probably  minimal  during  these  times.  It  seems  reasonable  to  select  a pan  constant 
of  1.0,  for  the  evaporation  from  the  pans  must  have  fairly  well  matched  that  from  the  pond,  be- 
cause the  effects  of  the  wind  shadow  and  slightly  elevated  pan  temperatures  tend  to  compensate 
for  each  other. 


Data  analysis 

Because  of  the  many  diverse  factors  affecting  evaporation  it  is  difficult  to  obtain  a precise  mathe- 
matical description  of  the  process  in  terms  that  will  allow  the  prediction  of  evaporation  under  spe- 
cific conditions.  (For  that  matter,  it  is  nearly  impossible  to  predict  with  any  great  confidence  what 
environmental  conditions  will  occur  at  any  time  in  the  future).  We  must  be  content  with  approxi- 
mations that  give  a range  of  values  for  evaporation  heat  loss  which  will  include  the  actual  value.  The 
data  spread  shown  in  the  evaporation  values  in  Table  II  for  the  various  investigators  in  this  field 
bears  this  out.  Even  the  well-controlled  laboratory  studies  of  Yen  and  Landvatter  (1970)  show  as 
much  spread  for  a water  temperature  range  of  0°  to  8°C  as  the  field  studies  at  Eielson  AFB  did  for 
a water  temperature  range  of  0°  to  1 6°C.  The  average  values  of  Ohtake  (1970)  show  a still  larger 
spread  for  a temperature  range  of  4°  to  1 5°C. 

Behlke  and  McDougall’s  values  (1973)  of  average  evaporation  compare  very  closely  to  those  of 
the  Eielson  experiments.  The  water  temperature  range  covered  by  their  experiments  was  2.7°  to 
1 1 .9°C;  their  average  value  of  evaporation  was  within  1 2%  of  both  Ohtake  and  the  Eielson  studies 
and  was  within  22%  of  Yen  and  Landvatter’s.  In  spite  of  the  spread  of  data,  it  is  interesting  to  note 
that  the  average  values  from  all  investigations  with  water  temperatures  below  16°C  fell  between 
3600  and  5698  g m'2  day"1 . This  is  a maximum  difference  of  only  45%.  Considering  the  differ- 
ences in  measurement  technique,  environmental  setting,  and  water  temperature  range,  the  average 
values  given  in  Table  II  are  considered  to  be  in  rather  good  agreement  and  probably  represent  the 
actual  magnitude  of  evaporation  under  unstable  conditions  ( Tw  > Ta)  rather  well. 

Comparison  of  the  measured  values  of  evaporative  heat  loss  to  the  various  formulas  proposed 
yields  some  interesting  results.  The  constants  o and  b in  Dalton’s  fundamental  evaporation  formula 
have  been  evaluated  for  the  Eielson  AFB  data.  To  do  this,  a computer  program  was  developed  to 
substitute  different  values  of  a and  b into  Dalton’s  equation  and,  using  the  wind  and  vapor  pressure 
data  for  each  measurement  at  Eielson,  to  calculate  the  theoretical  value  of  evaporation.  This  value 
was  then  compared  to  the  actual  measured  evaporation.  New  values  of  a and  b were  then  tried 
until  the  difference  between  the  calculated  and  the  measured  evaporation  was  minimal.  This 
was  done  for  each  of  the  78  data  points.  The  result  was  an  optimum  a and  b for  each  data  point. 
These  values  were  averaged  to  yield  an  equation  of  the  Dalton  form  which  gave  optimum  correla- 
tion with  the  Eielson  data.  The  optimum  value  for  a was  found  to  be  1 3.1 ; for  b it  was  found  to  be 
0.1 32.  Referring  to  Table  II  and  the  discussion  of  wind  during  ice  fog  conditions,  it  would  be  in- 
structive to  ipiore  wind  (which  was  absent  during  ice  fog)  and  calculate  the  value  of  a under  this 
restriction.  This  resulted  in  a value  of  a = 1 3.4.  The  standard  errors  of  these  two  equations  are 
within  1%  of  each  other,  confirming  that  the  effect  of  wind  during  ice  fog  is  minimal.  The  two  re- 
sulting equations  are 
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Eielton  Form  of  Dalton  Equation 
0I»(l3.l  + O.I32u)<*w-a.) 

Figure  5.  Evaporation  heat  loss  vs  Eielson  form  of  Dalton  equation 


Qt*  13.39  (•*-••) 

Figure  6.  Evaporation  vs  13.39  (ew-ej  (eq  17). 


Qv  = (1 3.1  +0.1 32w)(cw-ea) 


(W  m-J) 


= (13.4)  (ew~e  ) 


(Wr1). 


Behlke  and  McDougall  (197S),  who  made  a similar  correlation  for  their  data,  arrived  at  optimum 
values  fo  1 3.82  for  a and  zero  for  b.  It  is  interesting  to  note  how  closely  these  two  experiments 
agree.  Figures  S and  6 show  the  evaporation  data  plotted  against  the  Eielson  forms  of  Dalton’s 
equation  (eq  16and eq  17). 

Figures  7 and  8 show  the  evaporation  data  plotted  against  the  Kohler  equation  and  the  Russian 
Winter  equation.  The  Kohler  equation  undcrpredicts  the  measured  evaporation,  while  the  Russian 
Winter  equation  slightly  overpredic’s  it.  In  light  of  the  agreement  between  the  various  investiga- 
tors on  the  measured  values,  the  conclusion  is  that  currently  suggested  formulas  are  not  adequate 


Russian  Winter  Equation 


Figure  8.  Evaporation  heat  loss  vs  Russian  Winter  equation  (eq  7 1). 


for  describing  the  evaporation  under  the  highly  unstable  conditions  prevailing  during  ice  fog.  This 
is  probably  because  most  equations  fail  to  address  the  convective  activity  during  these  conditions, 
and  those  few  that  do  (primarily  the  Russian  Winter  equation)  use  empirical  values  obtained  from 
the  open  waters  of  rivers  and  lakes,  which  are  different  from  the  convective  activity  that  is  present 
over  the  warmer  water  of  a cooling  pond. 

The  Kohler  and  Penman  equations  (Kohler  1954,  Penman  1948)  contain  a wind  term  to  describe 
increased  evaporation  during  air  movement,  but  they  fail  to  account  for  the  fact  that,  although 
meteorological  wind  measurements  show  zero  wind  during  periods  of  ice  fog,  there  is  considerable 
air  movement  due  to  the  unstable  convective  conditions.  This  causes  evaporation  fluxes  that  are 
higher  than  the  equation  predicts.  The  Russian  Winter  equation,  on  the  other  hand,  is  much  closer 
to  the  actual  measured  conditions  and  even  overpredicts  the  measured  values.  This  suggests  that  a 
revised  torm  of  the  Russian  Winter  equation  type  might  be  formulated  that  will  give  a better  corre- 
lation than  either  the  Eielson  form  of  Dalton's  equation  (eq  16)  or  the  existing  Russian  Winter 
equation  (eq  11).  Figure  9 shows  the  data  plotted  against  an  equation  of  the  form: 
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Figure  10.  Evaporation  equations  vs  temperature  differences. 


Qe  = [4.84+(0.21)(rw-FJ)  (ew-ej  (Wm'2).  (18) 

This  equation  was  obtained  using  the  fitting  technique  discussed  above,  and  gives  the  optimum 
fit  to  the  data  from  the  Eielson  cooling  pond.  This  fit  considers  data  over  the  entire  temperature 
range  encountered  during  the  studies.  The  temperature  difference  (Tw-  Tt)  that  is  of  interest  dur- 
ing ice  fog  temperatures,  however,  is  that  of  AT  > 40°C,  since  the  pond  water  temperature  is  rela- 
tively stable  at  10°  to  1 2°C,  and  air  temperatures  of  at  least  -30°C  are  required  for  ice  fog.  This 
is  also  the  area  of  greatest  divergence  between  the  fitted  equations  and  the  Russian  Winter  equation 
(Fig.  10).  Restricting  attention  to  this  area  only  and  fitting  an  optimum  curve  (of  the  Dalton  form) 
to  the  data,  the  following  equation  is  obtained: 

<?E  = 14.21  (ew-e,)  (WnTs).  (19) 

This  should  represent  the  evaporation  from  a cooling  pond  during  ice  fog  conditions  rather  well, 
while  eq  18  would  give  a better  estimate  of  average  evaporation  for  general  conditions  during  win- 
ter months  in  interior  Alaska. 
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Figure  10  compares  all  of  the  equations  discussed  above.  It  could  be  argued  that  the  equation 
which  yields  the  maximum  values  (i.e.  the  Russian  Winter  equation,  eq  1 1)  should  be  used  for  calcu- 
lations of  evaporative  heat  loss.  This  can  be  supported  by  the  argument  that  the  extreme  instability 
supports  the  higher  evaporative  flux.  However,  the  greatest  instability  occurs  because  the  overlying  air 
is  at  its  coldest.  Under  these  conditions,  a greater  condensation  fallback  might  reduce  the  magnitude 
of  the  net  flux.  This  would  partially  explain  the  difference  in  the  constants  obtained  for  ice  fog 
conditions  (eq  18)  as  opposed  to  winter  conditions  in  general  (eq  11).  Equation  18  tends  to  over- 
predict the  measured  values  of  evaporation  at  AT  of  greater  than  40°C.  Therefore  eq  19  is  considered 
to  be  the  best  equation  for  severe  icv.  fog  conditions. 

SECTION  III.  RADIATION  DURING  ICE  FOG 

The  net  radiation  loss  during  ice  fog  from  an  open  water  area  such  as  a power  plant  cooling  pond 
is  dependent  on  several  variables.  The  emission  from  the  water  surface  is  dependent  on  the  tempera- 
ture and  emissivity  of  the  surface.  Although  surface  temperature  is  a difficult  parameter  to  measure 
with  high  accuracy  (Penman  1948),  it  can  be  measured  to  within  one  Celsius  degree  with  reasonable 
care.  The  emissivity  of  a water  surface  has  been  established  by  several  investigators,  and  the  value 
of  0.97  is  well  agreed  upon  (Raphael  1962,  p.  97).  Calculation  of  the  radiation  leaving  the  surface, 
therefore,  is  a relatively  easy  task. 

The  radiation  incoming  to  the  pond  can  be  divided  into  two  distinct  categories:  1 ) shortwave 
radiation  from  the  sun,  and  2)  longwave  radiation  from  the  atmosphere,  surrounding  geographic 
features,  or  man-made  structures. 

Shortwave  radiation  received  from  the  sun  is  a function  of  the  zenith  angle  of  the  sun,  rhe  tur- 
bidity of  the  atmosphere,  cloud  cover,  atmospheric  composition,  and  reflectivity  of  the  surface. 

Scott  (1964)  gives  the  maximum  daily  solar  radiation  arriving  at  the  earth’s  surface  under  clear  sky 
conditions  for  various  latitudes  in  the  northern  hemisphere.  On  a clear  December  22  at  the  latitude 
of  Fairbanks,  Alaska,  a daily  total  of  less  than  eight  langleys*  of  shortwave  radiation  is  received  at 
the  surface.  Koberg  (1964)  estimates  the  clear  sky  radiation  for  latitudes  up  to  50°N.  Extrapolat- 
ing his  computed  curves  yields  a value  of  less  than  10  langleys  of  radiation  reaching  the  surface  on 
21  December.  Behlke  and  McDougall  (1973)  calculate  a minimum  clear  sky  solar  radiation  reaching 
the  ground  at  the  latitude  of  Fairbanks  of  7 langleys  per  day.  Therefore,  the  magnitude  of  short- 
wave radiation  reaching  the  surface  during  December  should  lie  between  7 and  1 0 langleys. 

Clouds  and  ice  fog  will  reduce  this  value.  Eagleson  (1 970)  gives  an  extinction  factor  for  the 
effect  of  clouds  on  the  incoming  solar  radiation.  Table  III  lists  the  values  of  this  factor  for  various 
cloud  covers  and  cloud  heights.  During  ice  fog  the  cloud  height  is  zero  and  the  cloud  cover  is  total 
(10/10).  This  yields  an  extinction  factor  of  0.18.  Multiplying  this  factor  by  the  clear  sky  radiation 
received  at  the  surface  gives  a value  of  less  than  three  langleys  of  solar  radiation  received  per  day  by 
a horizontal  surface  during  ice  fog  at  this  time  of  the  year.  On  4 February  the  solar  radiation  arriv- 
ing at  the  surface  under  ice  fog  conditions  had  risen  to  slightly  less  than  10  langleys  per  day.  Com- 
paring this  to  the  average  measured  value  of  longwave  radiation  from  the  atmosphere  (430  langleys 
per  day,  Table  IV),  it  is  clear  that  the  incoming  solar  radiation  during  this  time  of  the  year  is  a 
negligible  portion  (~  2%)  of  the  radiation  heat  balance  of  an  open  water  surface  during  ice  fog. 

The  outgoing  longwave  radiation  emitted  by  any  body  can  be  calculated  from: 

<?„  = eoTo*  (20) 

* The  unit  of  solar  radiation  is  the  langley.  It  is  equal  to  1 g-cal  cm*2. 
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Table  III.  Cloud  attenuation  factors  [1-( 1-0. 18-0.024/ )N) 
for  shortwave  ra elation.  N = cloud  cover  in  tenths. 


Height  of  cloud  bm  z (ft) 


N 

0 

100 

1000 

4000 

8000 

16000 

0.1 

0.918 

0.918 

0.920 

0.928 

0.937 

0.956 

0.2 

0.836 

0.836 

0.841 

0.855 

0.874, 

0.913 

0.3 

0.754 

0.755 

0.761 

0.783 

0.812 

0.869 

0.4 

0.672 

0.673 

0.682 

0.710 

0.749 

0.826 

0.5 

0.590 

0.591 

0.602 

0.638 

0.686 

0.782 

0.6 

0.508 

0.509 

0.522 

0.566 

0.623 

0.738 

0.7 

0.426 

0.428 

0.443 

0.493 

0.560 

0.695 

0.8 

0.344 

0.346 

0.363 

C.421 

0.498 

0.651 

0.9 

0.262 

0.264 

0.284 

0.348 

0.435 

0.608 

1.0 

0.180 

0.182 

0.204 

0.276 

0.372 

0.564 

where  <?Q  = outgoing  longwave  radiation  (W  m"J ) 

6 = emissivity  of  the  radiating  body 
o = Stefan-Boltzmann’s  constant  (W  m'2  1C4) 

T0  = the  absolute  temperature  (0  K). 

The  problem  of  estimating  the  net  radiation  loss  during  ice  fog  conditions  is  therefore  reduced  to 
that  of  estimating  the  effective  back  radiation  from  the  atmosphere  and  subtracting  it  from  the 
outgoing  radiation  (eq  20).  Outgoing  radiation  is  fairly  straightforward;  however,  back  radiation 
is  not  so  simple,  and  several  investigators  have  proposed  empirical  equations  to  estimate  its  value. 

Brunt's  equation 

One  equation,  widely  used  for  atmospheric  radiation,  was  proposed  by  Brunt  (1944).  This  equa- 
tion uses  the  vapor  pressure  at  the  surface  to  estimate  atmospheric  radiation.  Approximating  condi- 
tions of  the  layers  of  air  at  various  heights  in  terms  of  the  conditions  at  the  surface  obviously  leaves 
much  to  be  desired,  but  can  be  justified  because  data  on  the  various  strata  of  air  above  the  ground 
are  seldom  available.  Anderson  (1952),  for  example,  found  that  twice  daily  radiosonde  balloon 
flights  were  not  adequate  to  accurately  measure  the  temperature  and  humidity  conditions  of  the 
various  strata  of  the  atmosphere.  This  equation  could  be  very  useful  if  it  can  be  shown  to  have 
sufficient  accuracy.  Brunt’s  equation  is 

Qa  = {a+bsfcUoTf)  (21) 

where  Qa  = atmospheric  radiation  (W  m ~2) 

ea  = vapor  pressure  of  the  air  at  2 m height  (mb) 
o = Stefan-Boltzmann’s  constant 
7,  = air  temperature  at  2 m (K) 
a and  b = empirically  derived  constants. 

Angstrom  equation 

Angstrom  (1920)  proposed  an  equation  similar  to  Brunt’s  but  with  slightly  different  form.  It  is, 
however,  somewhat  more  awkward  to  use.  Angstrom’s  equation  also  bases  the  prediction  of  in- 
coming atmospheric  radiation  on  the  vapor  pressure  of  the  air  near  tt>e  surface  - in  this  case,  at  4 m. 

Angstrom's  equation  is 

Qt  = [*-b  exp(-7ea)J(oft4)  (22) 
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Table  IV.  Longwave  radiation  data  from  ice  and  water  at  Eieison  AFB  cooling  pond. 


Net  longwave 


Date 

UtoWr  surface 
temp,  T 

ra 

Air  temp. 

7. 

ra 

Day  or 
night 

Fog 
cover * 

radiation 

outgoing 

(Wm-2) 

Average 

wind 

(mi-1) 

Over  water 

9 jam  73 

12.9 

-18.0 

D 

— 

89 

0 

12.9 

-18.0 

D 

102 

0 

18  Jan  73 

15.0 

-28.0 

N 

— 

159 

0 

20  Jan  73 

12.0 

-33.0 

D 

— 

106 

0 

21  Jan  73 

14.5 

-38.0 

N 

— 

102 

0 

22  Jan  73 

13.0 

-40.0 

D 

— 

82 

0 

23  Jan  73 

12.0 

-40.0 

N 

— 

122 

0 

27  Jan  73 

3.5 

-20.0 

N 

- 

63 

0 

5 Feb  73 

8.8 

-22.0 

D 

0.0 

112 

0 

7.5 

-14.0 

D 

0.0 

128 

0 

8.5 

-22.0 

N 

0.8 

114 

0 

6 Feb  73 

12.0 

-22.0 

D 

0.9 

106 

0.20 

12.0 

-18.0 

D 

0.0 

126 

0.20 

12.5 

-16.0 

N 

0.0 

128 

0.2C 

7 Feb  73 

13.5 

-20.0 

D 

1.0 

108 

0.05 

15.0 

-13.0 

D 

0.0 

142 

0.05 

8 Feb  73 

15.8 

-19.5 

O 

0.2 

152 

0.15 

16.6 

-16.0 

D 

0.2 

128 

0.15 

9 Feb  73 

16.7 

-23.0 

D 

0.5 

114 

0 

16.8 

-18.0  * 

N 

0.0 

162 

0 

10  Feb  73 

14.0 

-15.5 

D 

0.2 

158 

0 

11  Feb  73 

12.0 

-19.0 

D 

0.5 

159 

0.90 

13.8 

-16.5 

N 

0.0 

182 

0.90 

12  Feb  73 

12.3 

-18.5 

D 

0.0 

186 

3.0 

12.^ 

-23.5 

N 

0.0 

187 

3.0 

13  Feb  73 

13.0 

-26.5 

D 

0.2 

165 

1.50 

13.2 

-19.5 

D 

0.2 

184 

1.50 

14.0 

-17.5 

D 

0.5 

164 

1.50 

14  Feb  73 

12.7 

-32.0 

D 

0.9 

130 

0.25 

7 Mar  73 

7.5 

- 7.5 

D 

0.1 

65 

0.05 

8 Mar  73 

11.1 

-10.0 

D 

0.9 

77 

0.20 

11.1 

-10.5 

D 

0.9 

67 

0.20 

11.1 

- 9.5 

D 

0.9 

50 

0.20 

11.1 

-10.0 

D 

0.8 

76 

0.20 

9 Mar  73 

13.8 

- 9.5 

D 

0.5 

113 

0.30 

10  Mar  73 

14.4 

-13.0 

0 

0.1 

122 

1.60 

14.4 

-13.5 

D 

0.9f 

79 

1.60 

11  Mar  73 

14.5 

- 7.5 

D 

0.9f 

96 

0.15 

14.8 

- 7.0 

0 

0.0 

121 

0.15 

12  Mar  73 

14.5 

-15.0 

D 

0.4t 

137 

1.15 

14.5 

-15.0 

D 

0.3f 

109 

1.15 

14  Mar  73 

15.8 

-12.0 

D 

0.9 

90 

0.85 

16  Mar  73 

14.0 

-15.0 

D 

83 

1.30 

Table  IV  (cont'd).  Longwave  radiation  data  from  the  ice  and  water  at  Eielson  AFB  cooling  pond. 


Net  longwave 


Date 

Water  surface 
temp.,  Tw 

rc) 

Air  temp. 

J. 

rcj 

Day  or 
night 

Fog 

cover* 

radiation 
outgoing 
( Wm r*J 

Average 

wind 

(ms-'l 

19  Mar  73 

15.0 

- 5.6 

D 

144 

3.45 

21  Mar  73 

9.1 

- 9.0 

D 

0.0 

98 

0.10 

22  Mar  73 

6.5 

- 9.0 

D 

0.6 

103 

0.90 

Avg 

12.7 

-17.9 

0.4 

119 

0.60 

Over  ice 

27  Jan  73 

0.0 

-15.0 

N 

— 

66 

0.25 

28  Jan  73 

0.0 

-15.0 

N 

— 

66 

2.0 

-4.5 

-20.5 

D 

0.8 

60 

2.0 

-1.5 

- 1.0 

N 

0.9 

27 

2.0 

29  Jan  73 

-0.5 

- 1.0 

D 

0.9 

8 

3.0 

-1.0 

-10.5 

N 

0.9 

34 

3.0 

30  Jan  73 

-1.0 

-23.5 

D 

0.9 

49 

0.75 

-1.0 

-24.5 

N 

0.6 

42 

0.75 

31  Jan  73 

-1.5 

-25.0 

D 

0.9 

35 

0.15 

-1.5 

-25.5 

D 

0.0 

77 

0.15 

1 Feb  73 

-2.0 

-30.0 

D 

0.9 

52 

0 

-2.0 

-30.0 

N 

— 

34 

0 

2 Feb  73 

-2.5 

- 4.5 

D 

— 

73 

3.50 

-2.5 

- 9.0 

N 

— 

6 

3.50 

-2.5 

- 4.5 

D 

0.5f 

50 

3.50 

14  Feb  73 

— 

-20.0 

D 

0.0 

75 

0.25 

— 

-20.0 

D 

0.0 

79 

0.25 

— 

-19.0 

D 

0.0 

70 

0.25 

— 

-19.0 

D 

0.0 

68 

0.25 

— 

-17.0 

D 

0.0 

70 

0.25 

— 

-20.0 

D 

0.0 

70 

0.25 

— 

-20.0 

D 

0.0 

73 

0.25 

Avg 

-1.6 

-17.0 

0.4 

54 

1.20 

* Fog  cover  was  estimated  visually  by  the  observer, 
t Indicates  high  altitude  overcast. 


where  exp  = the  Napierian  base,  e = 2.72 
a,  b and  y = empirically  derived  constants. 

El  sasser's  equation 

Elsasser  (1942)  proposed  still  another  form  for  atmospheric  radiation  using  vapor  pressure  at 
4 m: 

Q,  = lp+b  log  et)  {oT*).  (23) 

The  above  equations  attempt  to  relate  the  atmospheric  radiation  to  local  vapor  pressure  only. 
While  these  make  convenient  equations,  they  are  obvious  oversimplifications  of  the  complex 
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meteorological  systems  governing  atmospheric  radiation.  A true  equation  should  be  related  to  the 
total  vapor  content  of  the  atmosphere,  the  carbon  dioxide  content,  the  ozone  content,  and  perhaps 
the  aerosol  content,  relationships  which  are  clearly  not  simple.  These  equations  can  only  be  looked 
upon  as  estimates  of  the  average  atmospheric  radiation,  and  evaluations  for  the  constants  involved 
in  each  equation  will  pertain  only  to  the  conditions  prevailing  for  that  set  of  data.  As  an  example, 
Anderson  (1952)  gives  values  for  the  various  constants  in  Brunt’s  and  Angstrom’s  equations  for  vari- 
ous localities.  Differences  as  great  as  135%  exist  between  different  studies  in  some  of  the  constants. 

Cloud  cover 

None  of  the  above  equations  considers  the  effect  of  cloud  cover.  To  compensate  for  this  effect, 
most  investigators  have  relied  upon  an  empirical  polynomial  equation  to  allow  back  radiation  to 
increase  as  cloud  cover  increases.  One  such  method  is  the  Karadag-Chumikova  factor  (Kondratyev 
1954).  This  technique  applies  a reduction  factor  to  the  net  loss  of  longwave  radiation  from  the  sur- 
face. From  observations  at  Karadag-Chumikova  in  southern  Russia,  it  was  found  that  a nonlinear 
polynomial  could  be  used  to  express  net  radiation  during  times  when  the  skies  are  cloudy.  Net 
radiation  was  then  expressed  as 

<?r  = Qn  (1-0.02/V-0.004/V2)  (24) 

where  QT  = net  radiation  under  cloudy  skies 
<?n  = net  radiation  under  clear  skies 
N = cloud  cover  in  tenths. 

A more  direct  approach  is  reported  by  Anderson  (1952).  He  found  in  studies  at  Lake  Hefner, 
Oklahoma,  that  a linear  equation  could  be  used  to  express  atmospheric  radiation  from  cloudy  skies. 
He  recognized  that  cloud  height  (in  addition  to  amount  of  cloud  cover)  must  be  taken  into  account. 
Anderson’s  cloudy  sky  equation  is 

Qac  = |[o+f>/Vexp{-0.058z)j  + [c-OAexp (-0.06z)]  (eg)|a7^  (25) 

where  <?ac  = atmospheric  radiation  from  cloudy  skies 
z = cloud  height  in  thousands  of  feet 
Ta  = air  temperature 

a,  b,  c,  D = empirically  derived  constants. 

This  equation  was  tested  for  cloud  heights  of  1600  ft  (500  m)  and  greater.  Anderson  states  that 
the  radiation  increases  as  cloud  height  decreases  and  probably  continues  to  increase  for  cloud  heights 
below  1600  ft.  He  contends  that  accuracies  of  5%  to  10%  are  possible  using  this  equation. 

Reflection 

The  amount  of  the  atmospheric  radiation  that  is  absorbed  by  the  surface  is  a function  of  the 
reflectivity  of  the  surface.  For  water  surfaces,  Raphael  (1962)  gives  a graph  of  reflection  as  a 
function  of  cloud  percentage  and  incident  angle.  For  angles  above  50°,  the  reflection  asymptotically 
approaches  0.05  for  overcast  skies.  Anderson  (1 952)  reports  that  studies  at  the  Physical  Standards 
Laboratory  of  the  U.S.  Geological  Survey  set  the  reflectivity  of  water  for  longwave  radiation  at  0.03. 
An  absorptivity  a of  0.97,  therefore,  seems  to  be  reasonable. 

The  net  radiation  heat  losr  under  clear  skies  during  ice  fog  temperatures  can  be  expressed  as 
<?n  = <?ou«  mo,  (26) 
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where  <?oul  = outgoing  longwave  radiation 

Qatmo,  = longwave  radiation  incoming  from  the  atmosphere 


Q„  = aewT*-(aaT*)(a+b)(a+by/Ta) 
= oew  T*  -( octT * )\a-b  exp  (-7<?a)] 
= oewT*-(aaT*)(a+b  logea) 


(Brunt  1944) 
(Angstrom  1920) 


(El  sasser  1942). 


The  choice  of  which  equation  to  use  for  atmospheric  back  radiation  is  not  critical,  since  as 

O 

Anderson  (1952)  points  out,  “The  difference  between  Brunt’s  and  Angstrom’s  equation  and  be- 
tween Brunt’s  and  Elsasser’s  equation  is  only  10%.  This  is  not  a significant  difference  since  the 
scatter  of  the  observations  is  considerably  greater.” 

Using  Brunt’s  equation  and  Kondratyev’s  Karadag-Chumikova  factor  for  cloudy  skies,  the  net 
radiation  Qr  can  be  expressed  as 

Q,  = °(ew^w(a+We7  )1  (1-0.02/V-0.004A/2)  (30) 

or  expressing  atmospheric  back  radiation  with  the  cloud  equation  proposed  by  Anderson  yields 
Qr  = aewT* -{{<7+b/Vexp(-0.0582)]  +[c-D/V  exp  (-0.06/)]  (ea)\aoT*  . (31) 

Comparison  of  the  above  equations  will  be  presented  later  in  this  section. 

Field  experimental  results 

At  the  Eielson  field  station  radiation  measurements  were  taken  several  times  daily  for  a period 
of  2'A  months,  from  January  to  the  middle  of  March.  Three  different  instruments  were  used:  a 
pyrgeometer,  a net  radiometer,  and  a precision  radiation  thermometer.  The  pyrgeometer  and  the 
net  radiometer  were  mounted  on  a cantilevered  beam  projecting  out  over  the  pond  surface  at  a 
height  of  0.5  m (Fig.  1 1 and  12). 

Radiation  values  were  measured  using  a Leeds  & Northrup  millivoltmeter  bridge  housed  in  a 
trailer  laboratory  at  the  site.  Calibration  of  the  instruments  was  performed  at  the  laboratory  of 
the  University  of  Alaska  Geophysical  Institute.  Measurements  were  made  at  frequent  intervals 
throughout  the  day  during  two  months  of  the  investigation  period.  The  longwave  pyrgeometer 
was  turned  toward  the  surface  and  then  toward  the  sky,  and  the  difference  between  the  two  values 
yielded  a "net  radiation”  which  was  compared  to  the  net  radiometer  readings.  Comparing  values 
from  the  two  instruments  showed  a correlation  coefficient  at  night  (when  all  incoming  solar  effects 
were  absent)  of  0.984.  As  spring  approached,  the  effect  of  incoming  solar  insolation  rapidly  in- 
creased each  day.*  A radiation  thermometer  was  used  to  measure  water  and  ice  surface  tempera- 
tures, and  to  calibrate  thermocouples  at  strategic  locations. 

The  average  value  for  outgoing  longwave  radiation  at  a water  temperature  of  10°C  averaged  10% 
below  the  theoretical  maximum  (Fig.  1 3).  This  is  probably  caused  by  the  attenuating  effect  of  fog 
between  the  instrument  and  the  water  surface,  and  the  field  of  view  ("shape  factor”)  of  the  instru- 
ment. Since  the  instrument  was  mounted  approximately  0.5  m above  the  water  surface,  the  water 
vapor  and  ice  fog  in  this  layer  attenuated  the  measured  outgoing  radiation  value. 


* The  radiation  pyrgeometer  used  at  Eielson  AFB  was  temperature  compensated  down  to  -25°C. 
Values  below  this  temperature  are  therefore  open  to  suspicion. 
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Figure  1 1.  Net  radiation  measurements  with  the  radiation  pyrgeometer. 


Figure  12.  Radiation  Instrumentation. 

Wexter  (1940)  measured  the  radiation  incoming  from  the  atmosphere  at  Fairbanks,  Alaska,  during 
the  winter  months.  Considering  only  values  below  0°C,  the  average  air  temperature  at  the  time  meas- 
urements were  made  by  Wexler  was  -1 5°C,  while  at  Eielson  AFB  the  average  was  -1 7°C.  The  cor- 
responding average  atmospheric  radiation  values  are 

Wexler  1 88  W m'2 

Eielson  214  W m*2 

Difference  10.8% 
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Figure  13.  Surface  temperature  t-s  outgoing  radiation  data. 


Foq  Cover 

Figure  14.  A verage  longwave  radiation  vs  fog  conditions. 

It  is  interesting  to  note  that  the  time  lapse  between  these  two  studies  was  34  years.  The  difference 
between  the  average  values  is  not  significant,  since  the  spread  in  both  Wexler’s  data  and  the  Eielson 
data  is  greater  than  10.8%.  Wexler  reports  values  from  1 32  to  279  W m'2,  while  the  Eielson  data 
contained  values  from  144  to  304  W m"2 . For  both  studies  this  represents  a spread  of  ±36%  about 
the  average.  Variations  of  as  much  as  69.7  W m'2  were  reported  by  Wexler  at  a specific  tempera- 
ture. Although  data  are  not  sufficient  to  be  conclusive,  the  overall  average  of  atmospheric  radiation 
appears  to  be  in  this  range. 

Analysis  of  experimental  results 

The  net  longwave  radiation  heat  loss  is  a function  of  the  atmospheric  conditions  prevailing  over 
the  pond.  The  equations  above  have  addressed  the  atmosphere  at  altitude;  however,  the  ice  fog 
plume  rising  from  the  pond  has  not  been  considered  by  these  equations.  The  back  radiation 
front  this  cloud  has  a direct  effect  on  the  net  radiation  heat  loss.  Figure  14  shows  the  net  radiation 
from  Table  IV  averaged  and  plotted  against  fog  cover.  This  illustrates  the  effect  of  fog  cover  on  the 
net  radiative  loss  of  heat  from  the  water  surface.  The  average  measured  radiation  loss  during  -30°C 
air  temperatures  (Table  IV)  is  37%  less  than  the  calculated  value  for  clear  skies,  using  the  corrected 
equation  discussed  below.  A corresponding  decrease  of  40%  is  found  when  net  measured  radiation 
is  compared  with  cloud  cover  (Fig.  14). 

Wien’s  Law  states  that  the  wavelength  of  maximum  intensity  of  radiation  emitted  by  a surface 
times  its  absolute  temperature  is  a constant  according  to  the  relation: 


TKux  ~ 0.2897  cm  K. 


If 


xy,: . ■ 


(32) 


The  average  temperature  of  the  water  surface  during  the  winter  months  was  12.6°C.  At  this 
temperature,  the  maximum  intensity  of  radiation  from  the  surface  occurs  at  a wavelength  of  10 
micrometers.  The  absorption  spectrum  of  water  vapor  has  a window  from  9 to  1 2 fim.  The  spec- 
trum of  the  radiation  is  wider  than  the  window,  so  that  a significant  portion  of  the  total  emission 
is  absorbed  by  water  vapor.  However,  approximately  30%  of  the  total  is  in  the  wavelength  range 
of  the  window  and  is  lost  into  space.  Geiger  (1965)  states:  “In  these  ranges,  streams  of  radiation 
can  escape  into  space.  ...No  matter  how  heavily  the  atmosphere  is  laden  with  water  vapor,  it  can 
afford  no  protection.”  Although  much  smaller  in  magnitude,  somewhat  similar  absorption  features 
exist  in  both  ice  and  liquid  water  absorption  spectrums  (Shaw  personal  communication  1974*, 
Kondratyev  1954).  (Water  has  much  more  intensive  absorption  bands  in  the  infrared.) 


Radiation  equations  for  ice  fog 


The  development  of  a precise  equation  for  the  net  radiative  heat  loss  from  an  open  water  surface 
is  extremely  difficult,  due  to  the  many  diverse  factors  affecting  the  various  components  of  the  radi- 
ation heat  balance  of  the  system.  Many  of  these  factors  cannot  be  measured  continuously;  others 
are  difficult  to  measure  accurately.  We  must  be  content  with  approximations  which  use  readily 
available  parameters  as  models  of  the  overall  system.  This,  of  course,  lacks  the  desirable  precision, 
but  in  this  manner  usable  average  values  and  their  bounds  can  be  calculated  for  use  as  design  para- 
meters. 
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Each  of  the  above  equations  for  atmospheric  radiation  contains  constants  which  must  be  deter- 
mined empirically.  Each  of  the  authors  has  determined  average  values  for  these  constants,  as  have 
other  users  as  well.  None  of  the  equations  address  the  situation  of  fog  immediately  over  the  water 
surface,  and  none  mention  ice  fog.  To  account  for  these  effects,  each  of  the  constants  was  evaluated 
using  data  from  the  Eielson  cooling  pond. 

To  evaluate  the  constants,  a computer  program  was  developed  which  substituted  various  values 
into  the  equations  and  compared  the  results  with  values  of  atmospheric  radiation  measured  at 
Eielson.  The  difference  between  the  calculated  and  measured  values  was  computed,  and  a search 
was  made  for  constants  that  would  reduce  the  difference  to  a minimum.  The  average  value  of  the 
constants  thus  derived  was  found  and  is  displayed  in  Table  V.  The  resulting  equation  was  again 
compared  to  the  measured  data,  and  a standard  error  calculated.! 

Although  the  differences  in  standard  error  are  not  great,  it  is  interesting  to  note  that  the  equa- 
tiors  originally  proposed  by  the  various  authors  generally  exhibited  considerably  greater  standard 
error  than  did  the  equations  resulting  from  the  Eielson  studies  (Table  V).  This  suggests  that  an 
increase  in  accuracy  of  8%  to  60%  can  be  obtained  when  calculating  atmospheric  radiation  during 
ice  fog  conditions  by  the  use  of  the  equations  correlated  to  the  Eielson  data. 

The  final  form  of  each  equation  is 

Revised  Brunt’s  eq  <?,  = (0.807+0.029  s/Tt  )oT*  (33) 

Revised  Angstrom’s  eq  (?,  = [0.658-0.201  exp(-0.375es))or4  (34) 


* Glenn  Shaw  of  the  Alaska  Geophysical  Institute, 
t Standard  error  is  defined  by  the  equation:  S.E. 


r(*cr*m«J2 
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Table  V.  Atmospheric  radiation  equations  (W  m*J ). 


i 


Std  error 

a b c d y Eg  no.  (96) 


Equation  constants  as  presented  by  original  author 


Brunt’s  eq 

0.60 

0.042  - 

21 

27.0 

Angstrom’s  eq 

0.806 

0.236  - - 0.115 

22 

24.3 

Elsasser’s  eq 

0.21 

0.22 

23 

73.8 

Anderson’s  eq 

0.74 

0.025  0.0049  0.00054 

24 

18.5 

Constants  resulting  from  Eielson  data 


Brunt's  eq 

0.807 

0.029  - 

33 

13.7 

Angstrom’s  eq 

0.658 

0.201  - - 0.375 

34 

14.6 

Elsasser’s  eq 

0.848 

0.039  - 

35 

14.1 

Anderson’s  eq 

0.814 

0.110  0.0054  0.00059 

36 

10.7 

Revised  Elsasser's  eq 
Revised  Anderson’s  eq 


Qa  = (0.848+0.039  log ea)oT* 

Qa  = [0.81 4+0.1 1/Vexp(-0.0584z)]+(eJ 

[0.0054-0.000594A/  exp(-0.06/)J  oT\ 


(35) 


(36) 


These  are  not  necessarily  the  best  or  most  representative  equations  that  could  be  found  to  repre- 
sent ice  fog  conditions.  However,  the  use  of  these  familiar  forms  allows  the  resulting  calculations 
to  be  compared  with  other  data  taken  in  different  localities  by  other  investigators.  The  standard 
error  of  10.7%  found  for  the  revised  form  of  the  Anderson  equation  (Table  V)  indicates  that  the 
results  of  calculations  using  this  equation  will  be  well  within  the  limits  of  the  measured  data  of 
both  this  and  other  studies.  Figures  1 5 through  18  compare  the  above  equations  with  the  data 
gathered  at  Eielson. 


Figure  15.  Atmospheric  radiation  data  vs  revised  Angstrom 
equation  (eq  28). 


29 


i 

-i 


k 


A 


Transmissivity  of  ice  fog 

It  was  noted  above  that  the  measured  outgoing  radiation  fell  10%  below  the  theoretical  value 
for  the  various  pond  temperatures  at  which  it  was  measured  (Fig.  1 3).  This  is  caused  by  two  factors. 
First,  the  measuring  probe  was  positioned  0.5  m above  the  water  surface,  in  this  position  it  re- 
ceived radiation  not  only  from  the  surface  of  the  pond,  but  also  from  the  much  colder  banks  of  the 
pond.  Therefore,  a radiation  “shape  factor"  should  be  applied  to  compensate  for  this  error.  In  addi- 
tion, the  0.5  m of  air  between  the  water  surface  and  the  probe  contained  water  vapor  and  ice  fog 
particles  which  absorbed  some  of  the  emitted  radiation  from  the  water  and  re-emitted  it  at  the 
temperature  of  the  vapor  and  ice  fog  crystals.  This  suggests  that  an  approximation  of  the  trans- 
missivity of  ice  fog  could  be  obtained  by  calculating  the  radiation  shape  factor  and  applying  it  to 
the  data;  the  difference  between  the  data  and  theoretical  value  would  then  represent  the  transmis- 
sivity of  the  ice  fog. 

The  radiation  probe  is  a hemispherical  dome  that  admits  radiation  from  all  angles  below  the  flat 
plate  upon  which  it  is  mounted,  f he  solid  angle  thus  represented  is  2n  steradians.  The  area  of  the 
pond  which  falls  into  the  “field  cf  view"  of  the  probe  represents  the  base  of  a cone  whose  vertex 
lies  at  the  center  of  the  probe.  The  solid  vertex  angle  of  this  quasi-conical  surface  represents  a good 
approximation  of  the  view  factor  of  the  radiometer.  This  can  be  calculated  by  integrating  the  plane 
vertex  angle  through  the  360°  field  of  view  of  the  probe.  Since  the  base  of  the  conical  surface  is 
difficult  to  express  as  a continuous  function,  the  technique  of  numerical  integration  can  be  used  to 
arrive  at  the  value  of  the  solid  angle.  The  resulting  value  of  the  shape  factor  is  0.956. 

A computer  program  was  developed  to  determine  the  difference  between  the  measured  values 
and  the  theoretical  value.  A correction  factor  for  the  theoretical  equation  was  found  which  gave  a 
minimum  standard  error  between  the  measured  data  and  the  formula.  This  value  was  found  to  be 
0.895,  and  is  a function  of  the  emissivity  e,  the  shape  factor  S,  and  the  transmissivity  r of  the  ice 
fog: 

-SBeBarBT  + efff7?  + -Swcwff7wr  ~ 0.895  oT*  (37) 

where  subscript  B refers  to  the  bank,  subscript  w refers  to  the  water  surface,  and  subscript  f refers 
to  the  ice  fog.  Rearranging  and  canceling  terms,  the  value  of  the  transmissivity  is  found  to  be 

0.895  T*  -ef  T? 

t = - . (38) 

All  of  these  terms  have  been  discussed  above  except  the  emissivity  of  the  banks  eB  and  the  emis- 
sivity of  the  ice  fog  ef.  The  emissivity  of  the  banks  can  be  assumed  to  be  that  of  snow,  e = 1 
(Eagleson  1970).  If  the  emissivity  of  the  fog  is  assumed  to  be  represented  by  any  one  of  the  re- 
vised atmospheric  radiation  equations  above  (eq  33,  34,  35  or  36)  an  order  of  magnitude  value  of 
0.8  is  representative.  These  equations  are  designed  to  measure  the  radiation  from  the  entire  atmos- 
phere. Since  the  thickness  of  fog  viewed  by  the  pyrgeometer  (300  m maximum)  is  a small  fraction 
of  the  thickness  of  the  atmosphere  (0.1  at  best),  an  emissivity  of  0.08  should  give  a reasonable  ap- 
proximation of  this  value.  An  approximate  value  of  the  transmissivity  is  then 


(0.044)  (243)4 +(0.97(0.956)  (283)4 


t = 0.9  (40 

constituting  an  absorption  of  10%  by  the  ice  fog  in  the  wavelength  range  of  emitted  radiation.* 


* This  assumes  that  the  density  or  the  ice  fog  is  constant,  and  can  be  represented  by  the  lowest  half 
meter. 
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Figure  19.  Corrected  radiation  data  vs  theoretical  outgoing  radi- 
ation (eq  29). 
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Figure  20.  Net  radiation  heat  loss  during  ice  fog  vs  tempera- 
ture difference  (Tw-  TJ. 

Figure  19  shows  the  data,  corrected  for  transmissivity  and  shape  factor,  plotted  against  the  for- 
mula of  eq  20. 

Figure  20  compares  the  four  equations  for  atmospheric  radiation  which  can  be  used  for  calculat- 
ing net  longwave  radiation  lost  from  the  water  surface. 


Radiation  model 

A model  of  the  radiation  balance  of  the  pond  during  ice  fog  conditions  may  be  constructed  in 
the  following  manner.  The  water  surface  will  be  considered  to  radiate  with  an  emissivity  of  0.97 
(Rout  = 0.97  oT*q  ).*  Some  of  this  outgoing  radiation  is  reflected  by  the  cloud  over  the  pond.  The 
ice  fog  cloud  radiates  at  the  rate  of  eoT*0.  This  is  a function  of  the  thickness  and  density  of  the 
fog  layer.  A portion  of  this  will  be  directed  back  to  the  pond.  Finally,  the  sky  above  the  pond  is 
radiating  longwave  energy,  some  of  which  is  directed  at  the  pond.  The  amount  reaching  the  sur- 
face is  a function  of  wavelength,  air  temperature,  transmissivity  of  the  ice  fog  cover,  and  shape 
factor  between  the  two  bodies.  If  an  approximate  effective  sky  temperature  7"JO  is  assumed,  then 
the  net  radiation  lost  Rntt  is 

* The  subscript  o will  be  used  to  indicate  absolute  temperatures. 
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^net  ~ ^outjoing-^reflected'^ic*  fog'^jky 

and  the  first  approximation  is 


(41) 


where  ef  = emissivity  of  fog 
ew  = emissivity  of  water 
Kt  = coefficient  of  reflection  of  ice  fog  cloud 
K(  = transmissivity  of  fog  layer 
5f(>  = shape  factor  between  the  fog  and  the  pond 
e,  = emissivity  of  atmosphere 
= absorptivity  of  water 
Twa  = water  surface  temperature 
7*,  = air  temperature. 

If  an  effective  sky  temperature  TK  is  assumed  to  be  a function  of  the  air  temperature  T^,  it  can 
be  seen  that  as  the  air  temperature  declines,  the  last  two  terms  of  the  equation  diminish,  thus  caus- 
ing  Rntt  to  increase.  As  fog  forms  over  the  pond,  however,  the  magnitude  of  the  radiation  from 
this  term  begins  to  increase  due  to  increasing  ef;  this  will  have  the  effect  of  attenuating  the  magni- 
tude of  net  radiation  lost,  and  an  inflection  point  in  the  curve  of  radiation  loss  vs  air  temperature 
will  result.  As  the  air  temperature  continues  to  decrease,  however,  the  decreasing  value  of  7^  will 
overshadow  the  increase  in  ef  (which  will  be  approaching  the  maximum  value),  and  a second  in- 
flection point  will  occur  after  which  the  net  radiation  lost  will  again  increase  with  decreasing  temp- 
erature. The  effective  sky  temperature  is  influenced  by  the  amount  of  water  vapor  in  the  upper 
layers  of  air  above  the  pond.  Pond  surface  temperature  also  varies.  The  degree  to  which  T%  is  a 
function  of  7~a  is  still  another  variable.  The  interaction  of  these  variables  causes  considerable 
scatter  (e.g.  see  Wexler  1940). 


Wind  effect 


The  effect  of  wind  (on  the  rare  occasions  when  it  is  present  during  ice  fog)  on  radiation  losses 
is  significant  since  it  can  greatly  reduce  or  completely  eliminate  the  ice  fog  cloud  and  its  back 
radiation.  When  the  wind  velocity  is  high,  the  vapor  laden  air  above  the  pond  is  swept  away  before 
it  forms  the  traditional  cloud  or  plume  over  the  pond.  Under  these  conditions  the  area  above  the 
pond  is  predominantly  air  from  the  surrounding  environment,  and  does  not  reflect  the  normal 
vapor  laden  conditions  present  when  wind  is  absent.  The  back  radiation  from  the  atmosphere 
under  these  conditions  will  be  considerably  lower  than  when  the  ice  fog  plume  rises  over  the  pond 
undisturbed.  It  is  evident  in  the  measurements  that  air  rrr  vement  played  a significant  effect  on 
the  net  radiative  loss.  Comparison  of  average  radiative  losses  on  windy  days  to  average  losses  on 
windless  days  in  the  same  temperature  range  shows  a 26%  higher  loss  on  windy  days  (Table  IV). 
Data  are  not  extensive,  but  it  does  appear  that  wind  (when  present)  plays  an  important  role  in 
radiative  as  well  as  other  heat  losses. 


Summary 


Ice  fog  plays  a significant  role  in  attenuating  the  radiative  heat  losses  from  an  open  water  sur- 
face. This  effect  is  more  pronounced  as  the  temperature  falls,  due  to  increased  back  radiation 
from  the  fog  layer.  Back  radiation  reaches  a maximum  and  then  declines  as  the  temperature  drops 
further.  Approximately  one-third  of  the  emission  from  a body  radiating  at  the  temperature  of  the 
cooling  pond  is  in  the  9 to  1 2 «m  wavelength  range  of  the  absorption  window  for  water  vapor.  Ab- 
sorption features  for  liquid  water  fogs  and  ice  particle  fogs  are  altered  but  still  exist  in  somewhat 
diminished  form. 
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Solar  radiation  can  be  ignored  for  December  and  January  in  interior  Alaska,  but  should  be  con- 
sidered during  daylight  hours  of  other  winter  months. 

The  outgoing  radiation  from  the  water  surface  at  1 1 °C  is  350  W m'1 . Outgoing  radiation  from 
the  ice  surface  at  an  average  -2°C  is  286  W m'2.  Using  the  revised  Anderson's  equation  (eq  36)  for 
the  back  radiation  from  the  ice  fog,  a value  of  162  W m*2  is  obtained  (for  -30°C  air). 

The  average  net  radiation  heat  loss  if  die  entire  pond  is  open  therefore  can  be  8.7x  1 06  W at 
-30°C.  This  is  48%  of  the  average  net  heat  input  to  the  pond  (see  Section  VI).  If  half  of  the  pond 
is  ice  covered,  then  24%  of  the  total  heat  is  dissipated  by  radiation  from  the  water  surface  and  1 5% 
from  the  ice  surface.  This  represents  a significant  portion  of  the  heat  budget  of  the  pond. 

SECTION  IV.  CONVECTION 

Heat  losses  by  convection  are  usually  the  largest  single  component  of  the  heat  budget  of  a body 
of  water  (Dingman  etal.  1967,  Rimsha  and  Donchenko  1958).  Convection  may  be  divided  into 
two  classifications:  laminar  and  turbulent.  However,  as  Gebhart  (1969,  p.  293)  points  out,  “Many 
of  the  natural  convection  processes  found  in  nature  are  of  such  a large  size  scale  and  of  such  long 
duration  that  the  detailed  transport  mechanism  is  largely  turbulent."  Visual  observations  of  open 
water  during  ice  fog  tend  to  support  this  point  of  view.  The  vapor  rising  and  condensing  over  the 
water  provides  an  excellent  visual  reference  of  the  local  eddies  and  overall  transport  of  the  vapor 
laden  air  near  the  surface,  which  is  seen  to  be  turbulent  with  laminar  flow  predominating  only  in 
the  first  few  centimeters  above  the  water.  Theoretical  studies  of  turbulent  heat  transfer  have  been 
made  by  many  researchers  (e.g.  Kays  1966,  Gebhart  1969,  Rotem  and  Ciaassen  1967).  However, 
only  Rotem  and  Ciaassen  have  addressed  the  subject  in  a manner  that  is  applicable  directly  to  free 
convection  heat  transfer  in  the  environment.  Their  approach,  although  rigorous  in  theory,  "found 
only  moderate  agreement  with  experimental  data."  Kays,  on  the  other  hand,  using  an  equally  rig- 
orous theoretical  approach,  demonstrated  good  agreement  between  theoretical  equations  and  lab- 
oratory experiments.  A question  arises  as  to  which  approach  to  use  - the  free  convection  approach 
of  Rotem  and  Ciaassen  (which  leaves  something  to  be  desired  in  verification),  or  the  forced  convec- 
tion approach  of  Kays  (which  correlates  well  with  conditions  that  may  be  far  removed  from  those 
encountered  in  the  environment). 

A third  alternative  is  available,  however,  and  that  is  to  rely  on  a purely  experimental  approach, 
as  did  Rimsha  and  Donchenko  in  deriving  an  empirical  formula  for  convective  heat  loss  from  data 
collected  throughout  the  Soviet  Union  in  the  winter  months. 

Other  empirical  approaches  have  attempted  to  determine  the  convective  heat  loss  by  direct 
measurement  of  the  volume  of  air  and  its  temperature  rise  while  passing  over  a pond  (Pruden  et  al. 
1954).  This  approach  has  the  advantage  of  yielding  values  that  are  not  dependent  on  other  vari- 
ables; however,  it  is  a difficult  and  time-consuming  process  to  use  and  does  not  lend  itself  to  ac- 
curate heat  loss  predictions. 

This  section  will  compare  the  various  approaches  and  discuss  their  correlation. 

Convective  heat  transfer  theory 

Consider  a plane  facing  upward  (e.g.  a water  surface  with  air  above).  The  minimum  convective 
heat  loss  (zero)  will  come  during  periods  of  zero  temperature  difference  AT  between  the  water 
and  the  air.  As  A T increases,  the  convective  heat  transfer  (?H  increases;  and  for  every  temperature 
difference,  there  is  a maximum  and  minimum  #H  depending  on  air  movement,  which  is  described 
by  the  convective  heat  transfer  coefficient  h of  the  air/water  boundary  layer. 


The  Prandtl  number  is  the  ratio  of  thermal  boundary  layer*  thickness  to  velocity  boundary 
layer  thickness.  When  these  two  are  of  nearly  equal  thickness  (Pi  = 1 ),  then  the  two  boundary 
layers  may  be  considered  equal.  Fluids  and  gases  with  Prandtl  numbers  approaching  1 can  be 
treated  as  though  their  boundary  layers  were  the  same.  For  these  cases,  the  momentum  and  energy 
equations  of  flow  can  be  assumed  identical;  this  greatly  simplifies  the  derivation  of  heat  transfer 
coefficients  and  equations. 

The  Prandtl  number  Pr  may  also  be  defined  as 

Pr  = cp/i/*  (42) 

where  cp  = specific  heat 

ti  = absolute  viscosity 
k = thermal  conductivity. 

For  air  at  temperatures  between  -35°  and  0°C,  this  turns  out  to  be  0.72. 

The  Grashoff  number  Gr  is  a dimensionless  number  that  has  been  found  to  characterize  free 
convective  flow,  much  the  same  as  the  Reynolds  number  characterizes  forced  flow.  The  Grashoff 
number  is  the  ratio  of  buoyancy  forces  to  viscous  forces  in  a free  convection  flow  system. 

The  Grashoff  number  Gr  may  also  be  defined  as 
_ gfSx'iK-TA 


where  g = acceleration  of  gravity 

0 = coefficient  of  thermal  expansion  = 1 / "F0 
x = distance  downstream  from  the  leading  edge  of  the  plane 
Tw-fg  = temperature  difference  across  the  boundary  layer 
v = kinematic  viscosity. 

The  Nusselt  number  Nu  is  the  ratio  of  the  length  along  the  surface  x to  the  boundary  layer 
thickness  8: 

Nu  = i = T = f(Pr‘ Gr)  (44) 

where  h = convective  heat  transfer  coefficient 
K = thermal  conductivity. 

It  is  the  function  of  the  product  of  Prandtl  number  and  Grashoff  number. 

The  convective  mechanism  prevailing  over  an  open  body  of  water  during  periods  when  much 
colder  air  covers  the  area  is  one  that  approximates  a single  large  convection  cell.  Air  near  the 
surface  of  the  water  is  heated  by  contact  with  the  water  and  by  latent  heat  of  evaporation  as  vapors 
condense.  The  heated  air  rises  due  to  the  buoyancy  forces  created  during  the  heating  process.  The 
rising  warm  air  is  replaced  by  cold  air  which  enters  the  area  around  the  periphery  of  the  pond,  lake, 
or  river.  This  gives  rise  to  a horizontal  flow  across  the  surface  of  the  water  (while  the  air  is  heating) 


* The  thermal  boundary  layer  is  defined  as  the  thin  layer  next  to  the  wall,  whose  thickness  extends 
to  the  interface  where  conductive  heat  transfer  is  equal  to  convective  heat  transfer. 


that  gradually  changes  into  a vertical  flow.*  Because  of  this  change  in  direction  of  flow  it  is  diffi- 
cult to  use  the  conventional  convection  parameters  which  involve  this  distance  x along  the  plate. 

It  is  customary  in  such  cases  to  define  a “characteristic  dimension  x"  to  be  used  for  computing  the 
Grashoff  and  Nusselt  numbers. 

Values  of  the  Grashoff  number  above  10*  are  usually  considered  to  be  turbulent  (in  environ- 
mental heat  transfer,  values  of  10s  are  sometimes  found  to  be  turbulent)  (Holman  1972).  A 
characteristic  length  for  turbulent  conditions  can  then  be  calculated  from  eq  43: 


" ymuK-rj  • 

For  air  at  -20°C,  over  water  at  10°C: 


■y  (9.8)  (1/253)  (30)  1 ' 

= 0.42  m. 

This  is  confirmed  by  visual  observations.  It  is  apparent  that  virtually  all  flow  under  ice  fog  condi- 
tions can  be  considered  turbulent. 

Rotem  and  Claassen  equation 

A theor  itical  approach  to  free  convective  heat  transfer  above  a horizontal  surface  was  presented 
by  Rotem  and  Claassen  (1968).  They  were  able  to  obtain  a similarity  solution  for  the  partial 
differential  equations  of  momentum,  continuity,  and  energy  by  considering  a semi-infinite  heated 
plate  facing  upward,  and  assuming  that  density  effects  were  negligible,  except  for  the  buoyancy 
terms,  for  a Prandtl  number  of  0.7.  Numerical  integration  was  used  to  solve  the  boundary  value 
equations,  and  the  results  were  compared  to  experimental  results  using  a “semi-focusing  color- 
schlieren  apparatus."  Their  experiments  indicated  that  a critical  Grashoff  number  for  turbulent 
transition  was  as  low  as  106.  They  concluded  that  the  horizontal  plate  is  much  less  stable  than 
the  more  widely  studied  vertical  and  inclined  plates. 

The  Rotem  and  Claassen  equation  is 

Nu="w^  (47) 

where  n is  the  power  of  the  temperature  difference  function  (6). 

That  is, 

e = Tx,-T%,  = cxn 
where  c = a constant 

x = the  distance  from  the  leading  edge  of  the  plate 
Tx  = the  surface  temperature  at  position  x. 

Figure  21  shows  a plot  of  Nusselt  numbers  vs  the  distance  x along  the  plate. 


* If  any  wind  is  present,  it  is  superimposed  upon  the  convective  flow,  and  changes  the  shape  of  the 
plume  accordingly. 
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Figure  21.  Nusselt  number  vs  distance  from  leading  edge  of  plate  (eq  47). 


Kays'  approach 

The  classic?!  approach  to  convective  heat  transfer  problems  is  presented  by  Kays  (1966).  Starting 
with  the  differential  equations  for  momentum,  energy,  and  continuity,  and  applying  them  to  the 
boundary  layer,  he  develops  generalized  equations  for  heat  transfer  by  convection.  For  the  external 
boundary  layer,  assuming  the  Reynolds  analogy  is  valid,  Kays  develops  the  following  equation  for 
convective  heat  transfer: 


St  = (48) 

[1+0.172  Re-0-1  (5  Pr+51  n (5  Pr+1)-14)[ 

where  St  = Stanton  number* 

Re  = Reynolds  number.! 

Kays  states  that  this  equation  is  a “little  awkward  to  use;  ...in  the  Prandtl  number  range  of  gases 
it  is  quite  well  approximated  by  the  following  more  simple  expression:” 

St  = 0 0295  Re'0  2 (49) 

Pr°« 


For  the  Prandtl  number  range  of  0.5  to  10,  the  various  assumptions  involved  in  deriving  these 
equations  have  been  shown  to  be  reasonable  by  experimentation.  Good  experimental  data  are 
available,  and  both  eq  48  and  49  are  shown  to  be  in  excellent  agreement  with  the  data  (Kays  1966). 

The  Stanton  number  can  be  written 


pcvum  k Re  Pr 

where  k = the  thermal  conductivity 

u_  = the  wind  velocity  of  the  unaffected  stream 
p = the  density 


(50) 


* The  Stanton  number  is  defined  as  Nu/Re  Pr. 
t Re  = umxtv. 


i 


■? 
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/»  = the  convective  heat  transfer  coefficient 
cp  = the  specific  heat. 

Combining  equations  49  and  50,  the  convective  heat  transfer  coefficient  can  be  written  as 

h = St(|) Re Pr  = 0.0295  (Re)0-8  (Pr)0,6(j).  (51) 

The  convective  heat  loss  is  then  calculated  from 

<?H  « t>{Tw-TJ  (52) 

or 

<?H  « 0.0295  (Re)°-8(Pr)0-8(i)  (Tw-Ta)  (53) 

where  Ta  ~ the  temperature  of  the  air  above  the  water 
Tw  = the  water  temperature  at  the  surface. 

Russian  approach 

Rimsha  and  Donchenko  of  the  Soviet  Union  have  reported  success  in  predicting  convective  heat 
losses  from  open  water  during  winter  months  in  the  Soviet  Union,  using  an  empirically  derived 
equation  (Rimsha  and  Donchenko  1958).  They  found  that  they  could  express  the  convective  heat 
loss  by  the  following  equation: 
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(?H  = [3.9  + 0.17  (rw  - Ta)  + 1.9 til  (Tw  - Ta)  (W  m‘2)  (54)  j 

where  u = the  wind  velocity  at  2 m (m  s'1 ).  j 

This  equation  has  two  advantages.  First,  it  is  supported  by  data  taken  under  winter  conditions,  i 

so  it  has  the  best  correlation  with  measurements  of  any  of  the  above  equations.  Secondly,  like  1 

Kays’  equation,  it  considers  the  effect  of  the  wind  and  is  therefore  more  general  than  the  free 


Figure  22.  Convective  heat  transfer  vs  temperature  difference 
(TW-TJ  at  zero  wlndspeed. 
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Figure  24.  Convective  heat  transfer  vs  temperature  difference 
(T  -7 '.)  at  a windspeed  of  2m  S'1. 


convection  equation  of  Rotem  and  Claassen.  Although  wind  conditions  are  probably  not  a 
significant  factor  during  ice  fog,  for  other  winter  conditions  the  effect  of  the  wind  must  be  con- 
sidered. 

Figures  22,  23,  and  24  show  the  three  equations  plotted  as  a function  of  the  temperature 
difference  (7"w  - Tt)  for  windsneeds  of  0, 1,  and  2 m s-1. 

Analysis 

The  three  equations  considered  differ  markedly  both  in  derivation  and  form,  yet  agreement 
between  them  for  the  lower  temperature  differences  is  good.  At  temperature  differences  of  as 
much  as  20  Celsius  degrees,  there  is  only  a 10%  difference  between  the  Rimsha  and  Donchenko 
and  the  Rotem  and  Claassen  equations  for  free  convection.  At  a 3 m s"1  windspeed,  the  Rimsha 
and  Donchenko  equation  and  the  Kays  equation  are  essentially  equal  for  temperature  differences 


V 


!200i 


Rimsha 

Oonchenko 


Temperature  Difference  (T*-T0) 

Figure  25.  Convective  heat  transfer  : •«  temperature  difference  at 
a windspeed  of  3 m S~l. 
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Figure  26.  Convective  heat  transfer  vs  temperature  difference  at 
a windspeed  of  5 m s~l. 


less  than  30  Celsius  degrees.  Since  the  Rotem  and  Claassen  equation  is  derived  for  free  convection 
and  contains  no  wind  term,  it  does  not  apply  to  the  cases  of  wind  greater  than  zero;  however,  it 
is  included  in  Figures  22,  23,  and  24  for  reference. 

The  empirical  approach  of  Rimsha  and  Donchenko  applies  throughout  the  entire  windspeed 
range.  Kays’  equation,  however,  does  not  begin  to  compare  favorably  with  the  others  except  for 
winds  over  2 m s-1 . This  is  to  be  expected,  since  Kays’  equation  was  derived  for  an  external 
boundary  layer  and  should  not  be  valid  for  conditions  of  free  convection  where  the  traditional 
boundary  layer  does  not  exist.  This  is  the  realm  that  Rotem  and  Claassen  addressed,  and  their 
equation  exhibits  fair  agreement  with  the  experimental  approach  of  Rimsha  and  Donchenko  when 
a value  of  n (in  eq  47)  = -0.4  is  used.  (The  isothermal  plate  situation  (i.e.  n = 0)  yields  values  of 
Qh  that  grossly  underpredict  Rimsha  and  Donchenko’s  equation.)  Above  2 m s-1 , however,  Kays’ 
equation  agrees  closely  with  the  Russian  empirical  formula  (Fig.  25  and  26). 


During  ice  fog  conditions,  the  wind  is  essentially  absent,  suggesting  the  use  of  the  Rotem  and 
Claassen  or  the  Rimsha  and  Donchenko  equations.  However,  since  a considerable  difference  exists 
between  these  two  equations  for  the  higher  temperature  differences  present  during  ice  fog,  another 
comparison  would  be  desirable  to  verify  one  equation  or  the  other.  The  Bowen  ratio  provides  such 
an  opportunity  for  comparison.  Bowen’s  ratio  can  be  written  (Bowen  1926): 

*B  = 2*  = (6.1)00 -*)phdx  (55) 

ew“e» 

where  <?H  = the  convective  heat  loss 
<?E  = the  evaporative  heat  loss 
P-  the  pressure  in  the  same  units  as  e 
rw  = the  water  surface  temperature 
Ta  = the  air  temperature 

ew  = the  saturation  water  vapor  pressure  at  the  water  surface  temperature  in  the  same 
units  as  P 

ea  = the  water  vapor  pressure  at  the  air  temperature.* 

A comparison  of  Bowen’s  ratio,  using  eq  55  for  conditions  over  open  water  at  10°C,  evaporation 
data  from  Section  II,  and  the  convective  heat  loss  predictions  of  equations  47  and  44,  is  shown  in 
Table  VI. 


Table  VI.  Bowen’s  ratio. 
Tw  = 10°C;  P = 1030  mb. 


Temperature 

difference 

T^-T, 


Vapor 

pressure 

difference 


Bowen's 

ratio 

R» 


Evaporative 

heat 


Bowen's  ratio 
Rimsha  and  Rotem  and 
Donchenko**  Claassen** 


0° 

10 

6.168 

40 

0.45 

0.45 

-10° 

20 

9.412 

1.375  1 

48 

1.04 

-20° 

30 

1.710  1 

55 

1.70 

1.40 

-30° 

40 

11.766 

2.140  1 

63 

2.60 

1.80 

-40° 

50 

71 

2.80 

* (?E  values  are  from  Table  I,  curve  averaged  for  QE  vs  A T W m'2. 
**  From  Figure  22. 


It  is  interesting  to  note  that  for  higher  values  of  temperature  difference,  Rotem  and  Claassen 
underpredict  the  Bowen  ratio,  while  Rimsha  and  Donchenko  overpredict.  This  would  suggest  that 
still  another  equation  might  be  found  that  would  compare  more  favorably  with  the  Bowen  ratio 
shown  in  Table  VI.  Figures  27  and  28  show  plots  of  two  other  equations  proposed  for  convective 
heat  transfer  in  the  environment.  Fishenden  and  Saunders'  equation  (Sutton  1953)  applies  only 
to  free  convection,  while  Devik’s  equation  (Devik  1932)  is  an  empirical  relation  that  was  derived 
from  measurements  in  Norway. 


* Saturation  vapor  pressure  can  be  used  with  very  little  error  at  ice  fog  temperatures. 


Figure  27.  Convective  heat  transfer  vs  temperature  difference  at 
zero  windspeed. 


Figure 


28.  Convective  heat  transfer  « temperature  difference  at 
windspeed  of  5 m s'*. 


Fishenden  and  Saunders’  equation  is 
Qh  = 2.5  (rw-ra)V4  (Wm-2). 
Devik’s  equation  is 


<?H  = 0.0213  Tm  \/u+Q3  (Tw-Ta) 


where 


+ 273  C. 


(WrrT2) 


(56) 


(57) 


(58) 
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Obviously  agreement  between  different  theories  and  studies  is  not  forthcoming,  since  neither 
of  these  equations  shows  a good  correlation  to  the  other  or  to  any  of  those  previously  discussed. 

Heat  budget  at  the  surface 

A heat  balance  at  the  water  surface  may  now  be  written: 

<?ta  = <?E+<?H+<?Rnef-  <59> 

Air  temperatures  Tt  of  -30°C  are  characteristic  of  average  ice  fog  temperatures.  The  surface 
temperature  of  the  cooling  pond  will  be  assumed  to  be  10°C.  The  average  evaporation  at  these 
temperatures  is  1 55  W m"2  (Fig.  10).  From  Section  III  the  average  radiation  heat  loss  for  these 
conditions  (using  the  revised  Anderson  equation  for  atmospheric  radiation)  was  found  to  be  188 
Wm-2  (Fig.  20). 

From  Figure  27,  Rotem  and  Claassen  predict  the  convective  heat  loss  to  be  300 Wm"2;  the 
percentages  are  given  in  Table  VII.  Other  values  in  Table  VII  are  computed  in  the  same  way.  Figure 
29  shows  the  components  of  the  heat  balance  during  ice  fog  plotted  against  temperature  differences. 


Table  VII.  Percentage  heat  budget  at  the  surface. 

Temp 


diff  Qe  Qh  Qr  Qj  hQ 


rc) 

(%) 

(96) 

(96) 

( WrrT 2; 

( Wm ~2  °C-'J 

10 

25.5 

25.5 

49 

157 

15.7 

20 

25 

38 

37 

319 

15.9 

30 

25.5 

43 

32 

490 

16.3 

40 

24 

47 

29 

643 

16.1 

Avg  16.0 


Figure  29.  Surface  heat  loss  components  during  ice  fog  vs  tem- 
perature difference  (TW-TJ  at  a water  temperature  of -10° C. 
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Since  all  three  heat  losses  are  a function  of  temperature  difference,  an  overall  heat  transfer 
coefficient  h0  can  now  be  defined  such  that 

<?T  = ho 

where  <?T  = total  heat  flux  (W  m*2).  It  is  interesting  to  note  that  the  value  of  ha  remains  nearly 
constant  through  the  range  of  temperature  differences  considered. 

This  overall  heat  transfer  coefficient  could  be  used  for  quick  calculations  of  total  heat  transfer 
from  an  open  water  surface  during  winter  conditions. 

Conclusions 

The  large  volume  of  experimental  data  reported  by  Rimsha  and  Donchenko  (1958)  from  many 
different  sites  in  the  Soviet  Union,  along  with  the  relatively  good  agreement  with  the  theoretical 
approach  of  Kays  at  5 m s_1  windspeeds  (Fig.  26),  lends  support  to  the  use  of  Rimsha  and 
Donchenko ’s  equation  for  convective  heat  transfer  during  windy  conditions.  In  addition  to  this, 
the  good  agreement  with  the  theoretical  approach  of  Rotem  and  Claassen  (1969)  during  free 
convection  conditions,  along  with  the  very  close  agreement  with  the  Bowen  ratio,  gives  support  to 
both  of  these  equations  for  free  convective  heat  transfer  for  temperature  differences  up  to  20  Celsius 
degrees.  When  ice  fog  exists,  high  temperature  differences  are  combined  with  very  light  wind.  This 
condition  needs  much  study  to  develop  conclusive  relations  for  convective  heat  transfer  during  these 
times. 

Convective  heat  loss  represents  a significant  portion  of  the  overall  heat  loss  from  an  open  water 
surface  during  ice  fog  conditions.  For  temperature  differences  of  30  Celsius  degrees  and  winds  of 
1 m s'1 , the  convective  heat  loss  according  to  the  Rimsha  and  Donchenko  equation  is  340  W m~2 
(Fig.  22).  If  the  entire  pond  were  open  water,  this  would  be  as  large  as  90%  of  the  heat  input  (as 
shown  in  Section  VI).  If  this  condition  should  exist,  the  high  rate  of  heat  loss  would  cause  an  in- 
crease in  ice  cover,  restoring  the  heat  balance  between  convection,  radiation,  and  evaporation. 

SECTION  V.  ICE  FOG  SUPPRESSION  TECHNIQUES 

Ice  fog  is  usually  confined  to  an  area  near  the  source,  and  boundaries  of  the  cloud  are  usually 
sharply  defined.  A short  distance  outside  the  cloud,  visibility  can  be  unlimited;  while  only  a few 
feet  inside  the  cloud,  it  may  be  less  than  1 5 m.  A typical  experience  of  pilots  during  winter  opera- 
tions in  Alaska  is  to  find  unlimited  visibility  all  along  their  route,  but  at  their  destination  to  find 
a mound  of  ice  fog  covering  the  airport,  making  landing  impossible.  (Often  visibility  down  through 
the  relatively  shallow  layer  of  fog  is  such  that  the  airport  and  surrounding  buildings  are  visible  from 
directly  above,  but  as  the  approach  for  landing  proceeds,  the  runway  is  observed  through  an  in- 
creasingly greater  thickness  of  ice  fog.  Consequently,  the  visibility  problem  becomes  such  that  the 
runway  disappears  during  the  final  stages  of  the  approach.)  This  leads  to  the  speculation  that  per- 
haps the  visibility  problem  could  be  eliminated  by  a means  of  controlling  the  location  and  path  of 
the  fog  cloud  so  that  the  runway  (or  other  critical  visibility  area)  is  outside  it.  At  the  Fielson  Air 
Force  Base  runway  an  abundance  of  clear  air  exists  outside  the  perimeters  of  the  base.  The  only  off- 
base  concentration  of  population  sufficient  to  cause  ice  fog  in  observable  quantities  is  in  the  Moose 
Creek  settlement,  to  the  northwest.  To  the  west  and  south  of  the  runway,  a supply  of  clear  air  is 
always  present. 
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Moving  clear  air  over  the  runway  to  exclude  the  ice  fog  offers  an  interesting  possibility:  a line 
of  blowers  or  fans  along  the  west  edge  of  the  runway  might  be  used  to  blow  the  ice  fog  away.  The 
fans,  properly  spaced  and  moving  air  at  the  correct  velocity  so  that  turbulence  and  crosswind  are 
not  a problem  for  landings,  should  eliminate  the  landing  visibility  problem.  Assuming  that  the 
length  of  the  runway  (14,000  ft  or  4267  m)  must  be  kept  clear,  that  a height  of  9.1  m (30  ft)  would 
be  sufficient,  and  that  a velocity  of  2 m s-1  (4  knots)  would  not  cause  a disruptive  crosswind  for 
aircraft  operations  but  would  be  strong  enough  to  overcome  the  density  gradient  flow  of  the  ice  fog 
(usually  in  the  order  of  1 m s'1),  the  power  required  to  move  the  necessary  air  would  be  approxi- 
mately 1 50,000  hp  (1 .1  x 108  W).  Ignoring  transmission  line  losses,  this  would  require  an  additional 
power  generation  of  over  1 10,000  kW  at  the  power  plant,  thus  compounding  the  ice  fog  problem  in 
that  area.  The  capital  equipment  cost,  both  in  air  handling  equipment  adjacent  to  the  runway  and 
additional  generating  capacity  at  the  power  plant,  would  thus  be  somewhat  discouraging;  however, 
the  biggest  problem  probably  would  be  one  of  stopping  the  creation  of  local  eddies  of  ice  fog  from 
flowing  behind  the  fans  into  the  negative  pressure  area  created  by  their  suction,  thereby  negating 
their  usefulness.  The  ducting  required  to  prevent  this  on  a scale  the  size  of  the  runways  makes  the 
capital  investment  prohibitive. 

A more  logical  and  efficient  approach  to  ice  fog  suppression  is  to  eliminate  the  contaminant  at 
its  source.  The  source  is  diverse;  therefore,  this  approach  also  has  its  limitations.  However,  since 
the  cooling  pond  is  a major  producer  of  ice  fog  (Benson  1965),  elimination  of  its  contribution 
would  be  a significant  step  in  the  overall  suppression  problem. 

Injection  wells 

In  eliminating  the  cooling  pond’s  use  during  the  winter  months,  water  could  be  drawn  from  the 
pond’s  make-up  wells  to  be  used  directly  for  cooling  and  then  disposed  of  in  an  injection  well 
downstream  from  the  intake  wells.  In  order  to  avoid  using  water  that  is  too  cold,  a recirculating 
loop  would  be  used  to  temper  the  incoming  well  water.  With  this  installation  in  operation,  the 
cooling  pond  could  then  be  isolated  and  allowed  to  freeze.  Since  no  heat  would  be  injected  into 
the  pond,  no  open  water  would  be  present  to  produce  ice  fog. 

Several  advantages  exist  for  this  proposal,  along  with  a few  disadvantages.  The  system  is  simple 
and  offers  a positive  elimination  of  ice  fog  from  the  pond.  It  would  require  very  little  manpower 
to  operate  and  the  basic  components  are  reliable.  However,  it  would  require  more  frequent  cleaning 
of  the  condenser  tubes,  due  to  heavier  mineral  deposits  from  using  raw  water  directly  from  the  wells, 
and  the  injection  wells  would  have  to  be  installed.  Further,  should  a pump  fail,  the  plant  would 
face  a cooling  crisis;  therefore,  backup  systems  would  have  to  be  provided.  The  present  cooling 
pond  could  serve  as  emergency  backup,  either  as  a water  source  or  as  a replacement  for  the 
injection  wells.  The  wells  would  have  to  comply  with  health  and  safety  regulations,  since  the 
potential  exists  for  injecting  contaminated  water  into  the  groundwater  aquifer  which  is  used  by 
residents  of  the  area  for  their  drinking  water.  Even  with  these  considerations,  however,  it  appears 
to  be  a more  practical  approach  than  the  fans  discussed  earlier  in  this  chapter. 

Latent  heat  storage 

Another  approach  to  eliminating  ice  fog  from  the  pond  is  to  maintain  an  ice  cover  over  the  pond 
while  the  pond  is  still  used  for  the  primary  cooling  source.  This  requires  that  the  heat  from  the 
power  plant  be  dissipated  by  conduction  through  the  ice  and  by  melting  of  the  stored  ice  reserves. 

The  advantage  of  this  plan  is  mainly  its  economy.  The  present  installation  could  be  used  with  only 
a minimal  modification,  and  all  equipment  necessary  for  the  modifications  would  be  available  from 
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surplus  U.S.  Air  Force  equipment.  The  water  available  for  cooling  would  be  of  better  quality 
(lower  in  mineral  content),  having  been  treated  by  passage  through  the  pond,  and  all  the  necessary 
backup  equipment  would  be  already  installed.  In  addition,  no  adverse  effects  to  the  groundwater 
aquifer  would  need  to  be  studied  before  the  system  was  installed  and  tested.  Section  VI  discusses 
this  technique  and  its  testing  and  performance  in  detail. 

Monomolecular  films 

Another  method  of  suppressing  the  ice  fog  from  the  pond  is  to  use  heavy  alcohols  for  suppression 
of  evaporation  from  the  surface.  This  has  been  successfully  used  in  a number  of  reservoirs  in  warmer 
climates  to  reduce  evaporation  losses  (Bean  and  Florey  1968,  Crow  et  al.  1969).  Success  has  been 
good,  by  all  reports,  with  evaporation  suppression  ranging  from  a low  of  6%  to  as  high  as  62%.  The 
heavy  alcohols  form  a monomolecular  film  over  the  water  surface.  Although  one  (mono)  molecular 
film  is  theoretically  all  that  is  necessary,  in  practice  it  is  found  that  several  layers  are  always  applied 
due  to  difficulties  in  application.  It  is  easier  to  apply  an  excess  of  the  relatively  inex|.ensive  chemi- 
cals than  to  risk  not  applying  enough  for  complete  coverage.  Several  methods  of  application  have 
been  tried,  including  dissolving  in  appropriate  solvents,  spraying  slurries  of  water  and  alcohol,  and 
hand  spreading  from  a boat. 

The  problem  of  heat  transfer  through  the  film  is  also  of  concern,  however,  since  the  function  of 
the  cooling  pond  is  to  dissipate  waste  heat  from  the  power  plant.  Since  there  is  no  ice  to  melt  with 
this  technique  and  little  evaporation,  heat  must  be  dissipated  another  way.  If  the  temperature  of 
the  pond  is  to  remain  relatively  constant,  conduction  and  storage  will  not  change.  Radiation,  how- 
ever, will  increase  significantly  since  the  water  surface,  no  longer  covered  by  ice,  will  radiate  at  a 
much  higher  flux  than  the  colder  ice  surface.  Convection  losses  will  also  increase,  for  the  same 
reason.  From  the  heat  budget  measurements  during  the  melting  experiments  (Table  VIII,  p.  64)  the 
average  heat  loss  to  melting  is  seen  to  be  6.4x  106  W.  Subtracting  the  radiation  losses  measured  over 
ice  cover  (Table  IV,  p.  23)  from  the  net  radiation  measured  over  open  water  gives  a net  increase  of 
radiation  for  the  uncovered  pond  of 

119-  54  = 65  (Wm'2). 

Estimated  net  convection  losses  (Fig.  23,  p.  39)  for  a 1 m s"1  wind  and  -30°C,  average  approxi- 
mately 520  W m~2 . Subtracting  the  losses  during  melting,  when  the  pond  was  ice  covered  (1  m s'1 
and  10°C),  gives  a net  increase  of  convection  losses  for  the  uncovered  pond  of 

520-90  = 430  (Wm"2). 

These  two  losses  add  up  to  495  W m~2 , It  therefore  requires  1 2,929  m2  of  pond  surface  to  replace 
the  melting  heat  loss  that  is  not  available  by  this  technique.  Slightly  less  than  30%  of  the  area  of  the 
pond  is  therefore  needed  to  dissipate  the  required  amount  of  heat  for  power  plant  operation  at 
these  air  temperatures.  The  remaining  70%  of  the  pond  is  available  as  reserve  cooling  capacity.  The 
use  of  monomolecular  films,  therefore,  appears  to  be  technically  feasible  for  this  application. 

To  test  the  effectiveness  of  the  evaporation  retardant  during  below-freezing  conditions,  two  evapo- 
ration pans  were  used  to  measure  evaporation  at  the  Eielson  cooling  pond.  After  establishing  an  evapo- 
ration base  line,  hexadecanol  was  added  to  one  of  the  pans.  Figure  30  shows  the  results.  Evapo- 
ration in  the  treated  pan  dropped  by  over  80%  during  the  first  period,  and  continued  to  remain  signifi- 
cantly lower  than  the  control  pan  until  the  test  was  suspended  1 20  hours  later.  The  average  suppression 
during  the  period  was  84%.  On  22  March  hexadecanol  was  spread  over  the  entire  pond  surface. 

Figure  31 , although  of  poor  photographic  quality,  shows  the  results.  A breeze  had  concentrated 
the  film  on  one  side  of  the  pond  (left  side  in  photograph).  The  boundary  of  the  film  can  be 
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Figure  30.  Evaporation  suppress  on  with  hexadecanol.  (Note:  hexa- 
decanoi  added  to  pond  surface  on  22  March  may  have  contaminated 
control  pan.) 


Figure  31.  Hexadecanol  film  boundaries  on  the  water  surface. 


seen  as  a sharp  line  between  the  surface  without  vapor  and  the  uncovered  water  surface  with  a 
normal  vapor  cloud  visibly  rising  from  the  surface.  This  also  points  out  one  of  the  major  problems 
to  be  resolved  in  using  these  films  for  evaporation  retardation,  i.e.  maintaining  the  integrity  of  the 
film.  Even  the  slightest  breeze  or  water  movement  will  move  the  film  to  one  side  of  the  pond, 
leaving  areas  of  open  water.  This  problem  was  not  apparent  in  the  pan,  partially  because  a large 
surplus  of  hexadecanol  was  available  to  maintain  the  film,  and  partially  because  of  the  size  of  the 
surface.  The  film  has  enough  internal  strength  to  keep  itself  intact  on  the  smaller  pan  surface.  This 
leads  to  the  speculation  that  the  problem  could  be  solved  with  a floating  grid  which  breaks  the  pond 
into  small  areas  on  which  the  film  can  maintain  a continuous  cover. 

One  other  problem  needs  to  be  considered:  that  of  bacterial  degradation  of  the  film.  Fiavo- 
bacterium  and  Pseudomonas,  among  others,  are  commonly  found  in  ponds  and  lakes,  and  reportedly 
thrive  on  the  film  (Chang  et  al.  1962).  As  the  film  is  applied,  the  population  of  these  bacteria 
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increases  sharply,  feeding  on  the  film  and  its  degradation  by-products.  As  the  population  of  these 
bacteria  increases,  the  film  life  decreases,  until  the  utility  of  the  film  is  destroyed.  Although  no 
observations  of  this  problem  were  made  at  Eielson,  two  factors  are  favorable  for  reducing  the  magni- 
tude of  the  problem  when  heavy  alcohols  are  used  for  ice  fog  suppression.  The  first  is  the  tempera- 
ture of  the  film.  During  ice  fog,  the  film  temperatures  are  significantly  below  those  found  during 
the  experiments  reported  in  the  literature.  Lower  temperature  should  do  much  to  suppress  the 
growth  of  the  bacteria.  Secondly,  the  pond  can  be  disinfected  by  either  adding  a bactericide  to  the 
hexadecanol,  or  by  adding  it  between  applications  of  the  alcohol  during  periods  without  ice  fog. 

The  technique  of  applying  hexadecanol  offers  several  advantages.  It  is  the  simplest  and  least 
expensive  of  any  approach  thus  far  investigated.  Chemicals  for  an  entire  ice  fog  season  would  not 
cost  more  than  a few  hundred  dollars.  Applications  can  be  done  by  one  person  in  an  hour  and  re- 
quire only  a small  boat  and  a spreading  device  such  as  a hand-held  lawn  seeder.  (Actually,  the 
chemicals,  which  come  in  a granular  powder  form,  can  be  spread  by  "hand  sowing"  if  the  seeder  is 
not  available.)  Application  need  be  repeated  only  when  the  film  is  depleted  during  the  ice  fog  period. 
No  manpower  is  required  between  applications  for  monitoring  or  operation  of  equipment. 

On  the  negative  side,  however,  several  items  need  to  be  resolved  and  a few  disadvantages  do  exist. 
First,  the  floating  grid  discussed  above  must  be  designed  and  built  and  the  design  problems  worked 
out.  Secondly,  the  bacterial  degradation  needs  to  be  examined  under  cold  weather  conditions,  and 
the  disinfecting  techniques  need  to  be  established. 

Other  heavy  alcohols  are  available  and  indications  are  that  some  of  them  may  be  even  more 
efficient.  It  appears  from  present  investigations  that  the  longer  the  alcohol’s  carbon  chain,  the  more 
effectively  it  suppresses  evaporation  (Dressier  1969).  In  this  respect,  hexadecanol  (Clb)  may  be 
surpassed  by  heptadecanol  (C,  7),  octadecanol  (Ci  *),  nonadecanol  (C»  #),  or  eicosanol  (C20)-  None 
of  these  other  chemicals  have  been  tried  under  ice  fog  conditions;  however,  during  tests  in  warmer 
(above-freezing)  temperatures,  they  are  reported  to  surpass  hexadecanol  by  a considerable  margin 
(Noe  and  Dressier  1969).  Evidence  also  exists  to  suggest  that  combinations  of  the  different  alcohols 
may  have  some  advantage  (Dressier  1969).  An  organized  testing  program  should  be  implemented  to 
determine  the  optimum  chemical  under  actual  conditions. 

Finally,  a disadvantage  to  the  technique  is  that  the  suppression  is  not  complete.  The  84%  sup- 
pression obtained  with  hexadecanol  may  be  optimized  to  an  even  higher  value  with  some  of  the 
other  chemicals  listed  above.  However,  some  evaporation  will  still  take  place  and  a percentage  of 
evaporation  will  still  produce  ice  fog.  This  should  not  be  of  too  great  a concern,  however,  since  the 
pond  will  be  reduced  to  the  category  of  a minor  ice  fog  producer  instead  of  the  largest  source  in  the 
area. 

The  technical  problems  and  unresolved  questions  concerning  this  technique  seem  to  be  much 
smaller  in  magnitude  than  any  of  the  other  techniques  investigated  thus  far.  It  certainly  warrants 
further  investigation  before  expensive  capital  investments  are  made  to  implement  any  of  the  other 
techniques  considered. 

Use  of  heated  water  from  the  cooling  pond 

Up  to  this  point  attention  has  been  focused  on  dissipating  heat  in  the  pond  in  a manner  that 
would  not  produce  ice  fog.  It  is  perhaps  prudent  to  recognize  the  possibility  of  practical  use  of  at 
least  a portion  of  this  heat  input.  Average  incoming  water  temperature  was  2 S°C.  This,  of  course, 
is  too  low  for  direct  use  in  heating  of  residences.  However,  the  magnitude  of  heat  input  (1 7.9x  10‘ 

W)  suggests  that  some  practical  use  might  be  found  (especially  in  light  of  the  current  energy  crisis). 
Such  possibilities  as  using  the  water  to  heat  hangar  floors,  parking  lots,  unheated  garages,  and  ware- 
houses should  be  analyzed  for  their  economic  implications.  Heat  pumps  and  turbines  using  Freon 


48 


cycles  should  be  investigated.  Perhaps  this  waste  heat  might  be  used  to  warm  the  solid  waste  land* 
fill  and  stimulate  the  biological  processes  there.  This  would  be  better  than  the  complete  waste  of 
heat  in  the  cooling  pond.  Economics,  of  course,  will  have  to  prevail.  However,  in  light  of  today’s 
concern  with  the  environment  and  the  awareness  of  dwindling  energy  supplies,  perhaps  new  costs 
wiii  make  some  of  these  proposals  economically  practical. 


SECTION  VI.  LATENT  HEAT  STORAGE  EXPERIMENTS 

Any  method  for  ice  fog  suppression  is  subject  to  two  overriding  considerations.  First,  the 
method  must  be  inexpensive  (i.e.  cost  less  than  the  cost  incurred  due  to  the  presence  of  the  ice 
fog).  Secondly,  the  method  must  be  applicable  to  use  in  the  extreme  cold  of  ice  fog  conditions. 
When  applied  to  power  plant  cooling  ponds  these  considerations  quickly  eliminate  a number  of 
large  scale  solutions  (such  as  enclosing  the  entire  cooling  pond).  One  technique  that  appeared  to 
warrant  investigation,  however,  was  to  cover  the  open  water  surface  with  a thick  layer  of  ice.  This 
not  only  greatly  suppresses  the  ice  fog,*  but  also  furnishes  a sink  of  latent  heat  that  can  absorb  the 
heat  output  of  the  power  plant  cooling  water.  Thus,  the  power  plant  can  continue  normal  operation 
without  any  negative  effect  from  the  suppression.  The  cost  of  the  operation  would  be  low  since  it 
would  use  normal  arctic  weather  to  provide  cooling,  and  the  only  energy  requirements  would  be 
the  pumping  required  to  enhance  the  freezing  process. 

Freezing  rates 

As  an  ice  sheet  grows  in  thickness,  all  the  heat  loss  necessary  for  freezing  must  be  transferred 
through  the  ice  sheet  to  the  cold  air  above.  A first  approximation  to  the  rate  of  growth  of  the 
sheet  can  be  made  by  equating  the  latent  heat  of  the  freezing  water  to  the  normal  conductivity  of 
the  ice  sheet.  This  yields  Stephan's  Law  for  ice  growth  (Pounder  1965): 

x = yfekD/Lp  = Cy/D  (60) 

where  x = the  thickness  of  the  ice  sheet 

k - the  thermal  conductivity  of  the  ice 
L - the  latent  heat  of  fusion  of  the  water 
p = density  of  water 

D = the  freezing  exposure  at  the  surface  of  the  sheet  in  degree-time 
C = y/2kJTp. 

The  ice  growth  rate  is  seen  to  be  proportional  to  the  square  root  of  the  freezing  exposure.  Thus, 
the  time  required  for  the  sheet  to  grow  increases  as  the  square  of  thickness.  When  building  ice  sheets 
for  suppression,  this  situation  is  intolerable;  an  ice  sheet  of  several  feet  must  be  formed  in  a matter 
of  days  and  then  grown  again  after  each  ice  fog  episode.  To  speed  up  the  process,  the  sheet  must  be 
thickened  by  freezing  on  the  top  where  each  new  layer  of  water  is  in  direct  contact  with  cold  air, 
thus  eliminating  the  insulating  effect  of  the  ice  (or  at  least  restricting  it  to  the  skim  ice  on  top  of  the 
layer  of  water).  To  accomplish  this  fast  thickening  process,  water  must  be  applied  to  the  surface 
of  the  ice  in  successive  layers  of  the  proper  thickness  to  freeze  at  an  optimum  rate. 


* The  vapor  pressure  over  ice  is  only  slightly  lower  than  over  water,  but  the  surface  temperature 
of  the  evaporating  ice  surface  is  several  tens-of-degrees  colder  than  that  of  open  water,  reducing  the 
evaporation  to  the  point  where  it  is  no  longer  visible. 


49 


Cooling  pond 

The  Eielson  power  plant  cooling  pond  (Fig.  32)  is  167  m wide  and  305  m long  with  a dike  15m 
wide  dividing  it  into  two  ponds  76  m wide  and  305  m long.  A channel  30  m wide  connects  the  two 
sides  at  the  east  end  of  the  pond.  During  normal  operation,  water  enters  the  receiving  pond  and 
proceeds  through  this  pond  around  the  east  end  of  the  center  dike  to  the  west  end  of  the  intake 
pond,  where  it  is  pumped  back  to  the  power  plant.  In  the  winter,  ice  forms  naturally  over  approxi- 
mately the  western  half  of  the  intake  pond,  leaving  the  remaining  portion  of  the  intake  pond  and 
all  of  the  receiving  pond  unfrozen  (3.48x10*  m2).  The  average  pond  depth  is  3 m,  giving  a volume 
of  1.04x10s  m3. 

An  8-in.  (20-cm)  diam  pipe  runs  along  each  side  of  the  receiving  pond  with  36  spray  nozzles 
spaced  every  12  m.  This  provides  a means  of  spraying  water  onto  the  surface  of  the  pond.  The 
header  was  supplied  with  water  from  beneath  the  ice  cover  by  a 6-in.  centrifugal  pump.  Valving 
was  arranged  to  permit  spraying  on  either  or  both  sides  of  the  pond.  The  header  on  each  side 
can  be  drained  separately  to  avoid  freeze-up.  A diversion  line  was  installed  from  the  west  end  of  the 
receiving  pond  to  the  west  end  of  the  intake  pond,  enabling  the  cooling  water  from  the  power  plant 
to  be  short-circuited  directly  to  the  intake  pond.  In  this  manner  the  receiving  pond  side  is  isolated 
from  hot  water  input  from  the  power  plant. 

Two  plastic  sheet  dams  were  installed  across  the  connecting  channel  at  the  east  end  of  the  dike 
(Fig.  33)  to  eliminate  density  gradient  circulation  of  warmer  water  from  the  intake  pond  to  the 
receiving  pond.  The  receiving  pond  was  therefore  free  of  all  heat  inputs,  so  that  freezing  could 
take  place  by  natural  means. 

As  soon  as  the  natural  ice  sheet  was  5 cm  thick,  spraying  would  commence  to  start  thickening 
the  ice  cover.  When  a layer  of  water  1 cm  deep  had  been  applied  to  one  side  of  the  pond,  the 
valves  would  be  switched  and  the  opposite  side  of  the  pond  sprayed.  At  the  end  of  spraying  for  the 
day  (usually  around  midnight),  a 3-cm  layer  of  water  would  be  built  up  on  each  side.  It  was  found 
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Figure  33.  Polyethylene  film  circulation  dam. 


that  3 cm  would  freeze  in  an  8-hour  period.  When  an  ice  fog  event  began  (i.e.  the  temperature 
dropped  below  -25°C),  all  spraying  would  cease.  The  diversion  of  the  water  to  the  intake  pond 
would  be  stopped  and  the  water  would  be  directed  back  into  the  receiving  pond.  This  would  then 
give  the  open  water  of  the  intake  pond  an  opportunity  to  freeze.  Ice  cover  on  this  pond  was  found 
to  form  in  approximately  25  hours  at  air  temperatures  of  -25°C.  When  the  intake  pond  was  frozen, 
the  ice  cover  of  the  entire  pond  was  established. 

Ice  fog  would  not  be  generated  from  this  source  until  the  warm  cooling  water  of  the  power  plant, 
injected  into  the  receiving  pond,  melted  enough  ice  to  form  a hole  in  the  ice  sheet  and  expose  open 
water.  The  time  required  for  open  water  to  appear  was  a function  of  the  air  temperature  and  thick- 
ness of  the  sheet  at  the  beginning  of  the  episode. 

Early  in  the  freezing  experiments,  it  became  apparent  that  the  system  had  some  inadequacies. 

The  pump  supplying  spray  water  was  not  producing  sufficient  pressure  to  spray  the  water  far 
enough  onto  the  ice  sheet.  The  flow  and  fre€Z6xharacteristics  of  the  water  at  the  working  temper- 
atures were  such  that  it  would  not  flow  more  than  1 2 m before  freezing.  The  result  was  that  a 
strip  along  each  side  was  thickened,  but  the  center  was  left  unflooded  to  thicken  by  natural  means. 

The  first  spray  nozzles  had  a 0.95-cm  diameter  opening.  It  was  found  immediately  that  the 
spray  from  these  nozzles  was  too  finely  divided.  At  -20°C  none  of  the  wate:  leaving  the  nozzle 
would  reach  the  ground  in  the  liquid  state.  Instead,  snow  was  created  which  built  up  in  large 
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Figure  34.  Snow-ice  buildup  around  0.95-cm  diameter  spray  nozzles. 

mounds  along  the  bank  (Fig.  34),  and  added  little  to  the  thickness  of  the  ice  cover.  After  experi- 
menting with  different  sizes  of  nozzles,  a 1.9-cm  diameter  was  found  to  be  optimum.*  Larger 
nozzles,  of  course,  had  lower  head  loss  and  a greater  flow;  however,  this  generated  a problem  of 
its  own.  The  larger  stream  impinging  on  the  sheet  of  ice  had  great  erosive  power,  which  was  partly 
due  to  mechanical  erosion  and  partly  due  to  melting  caused  by  the  4 °C  water.  The  stream  would 
erode  through  a 30-cm  ice  sheet  within  60  min.  As  soon  as  a hole  was  opened,  any  unfrozen  water  ( 

on  the  surface  would  drain  back  into  the  pond.  This  was  corrected  by  installing  a rotatable  nozzle  j 

(Fig.  35).  The  stream  was  then  changed  to  a new  position  every  1 5 to  30  min.  j 

Ice  building  ! 

The  rate  at  which  an  ice  sheet  can  be  thickened  is  a function  of  many  variables.  Outside  air 
temperature,  flood  water  temperature,  wind  conditions,  cloud  cover,  initial  ice  temperature,  and 
method  of  flood  water  application  all  combine  to  determine  the  rate  at  which  ice  will  form.  In 
previous  chapters  most  of  these  factors  have  been  discussed  in  some  detail.  The  only  variable  that 
has  not  been  discussed  (and  the  only  variable  that  is  usually  controllable)  is  the  method  of  applica- 
tion. 

Ideally,  the  water  should  be  applied  in  very  thin  layers  at  0°C.  As  latent  heat  is  removed, 
crystals  of  ice  quickly  start  to  form  and  the  layer  becomes  a slush  and  then  solidifies  into  a sheet  of 
ice  without  further  flow.  The  slush  has  very  poor  flow  characteristics,  and  if  flow  proceeds  during 
this  state,  a spongy  porous  ice  is  produced.  The  result  is  low  density  ice  with  lower  latent  heat 
capacity  and  relatively  thin  unflooded  areas  between  flooded  sections  — an  undesirable  configuration. 

On  the  other  hand,  when  the  water  is  applied  well  above  freezing,  the  heat  content  of  the  water,  as 

* An  additional  problem  with  the  smaller  nozzles  was  their  tendency  to  freeze,  stopping  the  flow 
through  that  portion  of  the  pipe.  When  this  happened  at  the  end  of  the  header,  the  water  in  that 
portion  of  the  header  pipe  would  freeze  and  it  would  be  lost.  The  larger  nozzles  showed  no  tend- 
ency to  freeze  even  when  used  at  -50°C. 
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Figure  35.  Rotatable  nozzle. 


well  as  the  latent  heat  of  fusion,  must  be  removed.  A top  skim  of  ice  forms  first,  and  thereafter 
the  rate  of  freezing  is  subject  to  the  limitations  of  eq  60.  (Water  nearly  always  continues  to  flow 
out  from  under  this  top  skim  ice,  leaving  an  insulating  air  film  between  it  and  the  ice.  This  is  also 
an  undesirable  condition.) 

Nonuniform  thickness  is  undesirable  since  it  allows  subsequent  melting  to  create  areas  of  open 
water  and  form  ice  fog  before  the  bulk  of  the  latent  heat  stored  in  the  thicker  portions  of  the  sheet 
is  used.  Therefore,  a compromise  must  be  made  in  the  freezing  technique.  This  compromise  will 
vary  with  each  change  in  the  weather  conditions,  and  is  therefore  difficult  to  optimize.  Several 
methods  of  application  were  tried  and  each  one  will  be  discussed  separately.  However,  all  methods 
are  subject  to  the  maximum  amount  of  heat  that  the  system  as  a whole  can  remove. 

Ice  growth  rate  — maximum  system  limitation 

The  active  air  mass  that  covers  the  pond  is  approximately  7.5  m deep.  The  intake  pond  area  is 
2.3xl04  mJ,  giving  an  effective  air  volume  of  1.7x10s  m3.*  An  ice  growth  rate  of  2.54  cm  hour'1 
over  the  entire  receiving  pond  surface  requires  a pumping  rate  F of 

F = (2.3x104)(2.54x10-2)i.  = 9.7m3  min-'.f 

bU 

The  heat  Q that  must  be  removed  to  freeze  an  area  of  2.3  x104  m7  of  4°C  water  2.54  cm  deep  is 
Q = (mass  of  water) ((temperature  drop)(specific  heat)  + (latent  heat  to  be  removed)] 


* This  is  conservative  (on  the  high  side)  since  the  air  mass  affected  cannot  be  much  greater  than  the 
highest  arch  of  the  spray.  This  is  generally  below  6 m. 

f The  pump  installed  provided  a maximum  of  1 .5  m2  min'1  with  a discharge  head  of  5.17x1 02 
N m~2.  However,  it  could  be  set  up  to  provide  3.78  m2  min-1  if  only  2.76x10s  N m“2  were  re- 
quired. 
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Q = (2.3x10*  ml)(2.54x10'J  m hr_l)(1  xIO*  g m-3) 

((4°C)(1  cal  g-‘  °C~1)  + (80  cal  g~* )]  = 4.9x10‘°  (cal  hr1)  = 56.9x10*  W. 

If  we  assume  that  the  effective  air  mass  passing  over  the  pond  can  be  raised  5.5°C  (the  average 
recorded  temperature  rise  is  4°C,  upstream  to  downstream),  then  the  mass  of  air  necessary  to 
dissipate  this  heat  is 

_ _ <?  _ 4.9 xIO10 

cpAT  " (0.24)(5.5) 

= 3.7xl010  (ghr~*) 

= 2.5 xIO7  (m3  hr-1). 

The  effective  air  volume  is  1 ,7x  10s  m3 ; therefore,  1 82  air  changes  per  hour  are  needed.  If  wind 
blows  across  the  shortest  dimension  of  the  pond  (76  m),  then  a wind  velocity  of  (182)  (76)  = 
13,832  m hr-1  = 3.8  m s'1  would  be  required  to  form  2.54  cm  of  ice  per  hour. 

Evaporation  heat  loss  - maximum 

Heat  loss  to  evaporation  will  be  limited  by  the  temperature  and  incoming  humidity. 

Temp  Weight  of  water  in  saturated  air 

-20°C  2.6x10~4  kg/kg  of  dry  air 

-10°C  4.6x10“4  kg/kg  of  dry  air 

0°C  7.9  x 1 0"4  kg/kg  of  dry  air. 

If  we  assume  that  air  enters  the  pond  area  with  relative  humidity  equal  to  zero  and  leaves 
saturated,  and  the  effective  air  mass  flow  is 

Ma  = 4.04 xIO7  kg  hr'1 

then  the  maximum  possible  mass  and  heat  lost  to  evaporation  are 


Mass 

Heat  loss 

Temp 

(kg  hr') 

(W) 

~20°C 

10,504 

6.6  xIO6 

-10°C 

18,584 

11.6x106 

0°C 

31,916 

19.7x10*  . 

This  heat  loss  will  freeze  the  following  amounts  of  ice: 
Temp 

-20°C 
-10°C 
0°C 


kg  hr*1  cm  hr~‘ 


71,250  0.279 

125,000  0.543 

212,500  0.924. 


54 


Realistic  growth  rate 

Using  the  average  meteorological  values  measured  at  Eielson  of  A 7"  * 4.0°C,  average  wind  = 

1.7  m s“*  at  6 m above  ground,  and  an  effective  air  mass  of  1 .5  x 10s  m3;  and  assuming  that  ail 
heat  transfer  is  done  during  the  flight  of  the  droplet  before  impact  and  that  the  convective  heat 
transfer  efficiency  t)  is  ij  = 0.6,  the  rate  of  freeze  for  the  system  can  be  estimated.  The  convective 
heat  transfer  capacity  is 

Q„  = (T}(mcpA7)l  = (0.6)<1.5x10s)(1.7)(3600){0.24)(1400)(4)  = 11.3x10*  W. 

The  realistic  evaporative  heat  transfer  is  probably  about  30%  of  the  maximum  values  calculated 
above.  A maximum  realistic  value  would  then  be 

<?E  * S.9x106  W. 

Radiative  heat  loss  from  Section  III  would  be  84  W m~2.  Estimation  of  the  average  spray  stream 
area  at  1000  times  that  of  the  emerging  jet  approaches  realism,  since  the  streams  were  compact  and 
did  not  break  greatly  over  most  of  their  flight.  The  area  for  18  streams,  1.9  cm  in  diameter  and 
18  m long,  is 

As  = (18)(1.9xi0-2)(ir)(18)(1000)~20,000  m2. 

This  gives  a radiative  heat  loss  of 

<?R  = (84)(20,000)  = 1 ,68x10s  W 
and  the  total  heat  transfer  (?T  will  be 

QT  = <?H+QE  + <?R  = 11.3x106  + 5.9x106  + 1.68x106  = 18.8xl06  W. 

This  will  freeze  the  following  volume  of  4°C  water: 

V * = 19x105  kghr-i  = 190m3  hr' 

or  for  a uniform  spread  over  the  pond,  this  is  a growth  rate  G of 


2.3  xIO4 


100  = 0.83  cm  hr-1. 


This,  of  course,  ignores  the  considerations  of  convective  heat  transfer  from  the  droplet  to  the 
air,  as  well  as  the  heat  transfer  within  the  droplets.  It  is  assumed  that  these  are  not  the  limiting 
factors. 

It  is  easily  seen  that  under  the  above  conditions,  1 cm  hr'1  of  freezing  is  optimistic,  and  0.75 
cm  hr"'  is  probably  more  realistically  close  to  the  upper  limits  of  the  system. 

Peyton  and  Johnson  (1967),  in  their  Kotzebue,  Alaska,  trials,  were  only  able  to  sustain  a growth 
rate  of  8 in.  in  24  hours,  or  0.76  cm  hr-* . Johnson  feels  that  this  is  the  upper  limit  of  the  growth 
rate  when  spray  techniques  are  used.  (In  addition,  the  area  in  which  their  experiments  were  con- 
ducted is  subject  to  much  higher  wind*,  ban  are  normal  for  the  Eielson  area.) 
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Figure  36.  Intelliglant. 


Ice  building  techniques 


Spraying.  The  use  of  a spray  to  apply  water  to  the  surface  of  the  ice  sheet  has  the  advantage  of 
cooling  the  stream  of  water  before  it  strikes  the  ice.  The  amount  of  this  convective  cooling  can 
range  from  complete  freezing  of  the  stream,  so  that  snow  is  produced,  to  merely  cooling  the  water 
without  any  latent  heat  extraction.  The  degree  of  cooling  can  be  controlled  by  nozzle  configura- 
tion. Variable  orifice  nozzles  have  long  been  available  that  adjust  from  a thin  compact  stream  to  a 
very  finely  divided  mist.  The  cost  of  such  nozzles  is,  however,  prohibitively  high  for  a multi-station 
header.  Thus,  for  the  cooling  pond  experiment,  a fixed  orifice  nozzle  was  chosen.  The  stream 
could  still  be  varied  by  changing  nozzles  to  fit  the  atmospheric  conditions.  However,  at  ice  fog 
temperatures  the  chore  of  changing  36  nozzles  was  not  looked  upon  with  enthusiasm,  and  after  a 
few  changes,  a good  compromise  nozzle  was  instilled  and  rarely  changed  thereafter.  A 1 .9-cm- 
diameter  fixed  orifice  nozzle  was  found  to  give  good  compromise  between  distance  and  stream 
breakup  for  the  pressure  produced.  The  nozzle  configuration  consisted  of  a pipe  plug  with  a hole 
drilled  in  the  center.  The  inside  of  the  plug  was  naturally  curved  so  that  some  streamlining  of  the 
flow  occurred.  These  nozzles  were  tested  against  a group  of  commercial  nozzles  and  found  to  give 
comparable  performance  at  a fraction  of  the  cost.  (Fire  hose  nozzles  were  also  tested  and  found  to 
spray  approximately  1 Zi  times  the  distance  of  any  other  nozzles  tested,  but  at  a cost  of  several 
times  that  of  the  inexpensive  model  chosen.) 


The  headers  along  the  two  sides  of  the  receiving  pond  were  275  m long  and  contained  36  nozzles 
(1 8 on  each  side).  Two  45°  elbows  were  installed  at  each  nozzle  port.  The  elbows  provided  a 
means  for  rotating  the  nozzles  about  their  vertical  axes  to  give  wider  coverage  and  alleviate  the 
erosir-o  problem  at  impact  (Fig.  35). 

A second  spray  technique  involved  the  use  of  an  Intelligiant  (Stang  Hydronics,  Inc.)  shown  in 
Figure  36.  This  machine  is  basically  a targe  spray  nozzle,  with  its  direction  controlled  by  two 
hydraulic  cylinders.  One  cylinder  controls  rotation  about  the  vertical  axis  (azimuth)  and  the  other 
one  moves  the  nozzle  about  the  horizontal  axis  (elevation).  The  rate  of  cylinder  movement  and 
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the  throw  of  each  cylinder  are  adjustable,  giving  a very  versatile  spray  program.  The  cylinders  are 
driven  by  water  from  the  pump  which  supplies  the  nozzle,  so  that  a water  pump  is  the  only  power 
source  needed.  The  nozzle  used  is  a firehose  nozzle,  giving  the  device  good  range  and  excellent 
coverage.  A radius  of  30  m and  a sweep  of  90°  were  obtainable  with  this  unit,  providing  the  best 
coverage  of  any  spray  technique  tested  (730  m2  per  setting).  Spray  flow  rate  was  0.46  m3  min'*. 

Flooding.  Several  techniques  were  explored  for  flooding  the  surface  of  the  ice  sheet.  All 
suffered  from  the  inability  to  spread  the  water  evenly  over  large  areas  of  the  surface  of  the  ice  sheet, 
and  the  cooling  realized  from  the  spraying  technique  was  absent.  However,  in  some  cases  0°C  water 
could  be  drawn  from  just  beneath  the  ice,  tending  to  compensate  for  the  loss  of  convective  cooling 
during  spraying.  A short  description  of  various  flooding  experiments  follows. 

1 ) A submersible  pump  was  lowered  into  the  water  through  a small  hole  in  the  ice.  The  dis- 
charge was  directed  through  a T-pipe  onto  the  ice  surface.  The  requirement  for  electrical  power 
and  the  long  power  lead  reduced  the  versatility  of  the  pump  considerably.  The  discharge  was  too 
concentrated,  and  therefore,  erosion  was  also  a problem.  Nor  was  the  pressure  head  great  enough 
to  permit  a spray  discharge.  Pump  capacity  was  approximately  0.1 2 m3  min”1  (32  gal.  min”1 ) for  a 
one-half  horsepower  pump.  A larger  model  of  this  pump  in  conjunction  with  a portable  generator 
might  have  performed  better. 

2)  A common  ice  auger  can  be  used  as  a pump.  Its  pumping  rate  is  approximately  3 m3min”' 
(800  gal.  min"' ).  It  has  the  advantage  of  being  portable  so  that  flooding  can  be  provided  wherever 
more  ice  is  needed.  An  additional  advantage  is  that  only  one  person  is  needed  to  operate  it.  It  is 
fast,  since  it  drills  its  own  hole  and  then  pumps  out  of  it.  A disadvantage  is  that  it  sprays  the 
operator  to  the  knees,  so  he  must  wear  insulated  hip  boots.  Another  problem  is  that  the  weight  of 
the  water  on  top  of  the  ice  deflects  the  ice  sheet  downward,  so  that  a conical  basin  of  water  is 
created  on  the  surface  with  the  auger  hole  in  its  center.  Much  of  the  water  therefore  drains  back 
into  the  pond.  This  disadvantage  can  be  countered  by  a simple  sheet  metal  drain-back  shield  in  the 
form  of  a cylinder  20  to  25  cm  high.  This  shield  is  frozen  in  place  before  flooding  and  blocks  the 
drain-back  of  flood  waters  into  the  hole. 

If  the  thickness  of  the  ice  sheet  is  less  than  the  distance  between  convolutes  of  the  auger 
(approximately  20  cm),  the  auger  will  pull  into  the  hole  by  wedging  the  sheet  between  two  of  the 
convolutes  and  driving  itself  into  the  hole  as  the  auger  turns.  In  this  situation  a sleeve  of  adequate 
length  must  be  inserted  into  the  hole  to  eliminate  the  problem. 

3)  Outboard  motors  make  excellent  pumps,  and  will  flood  an  ice  surface  with  up  to  7.56  m3 
min”1  (2000  gal.  min”1 ).  They  have  the  same  problem  with  drain-back  as  the  auger,  and  they 
require  a larger  hole.  A hole  60  cm  in  diameter  will  suffice.  It  can  be  easily  cut  with  a chain  saw 
(Fig.  37).  The  pool  of  flood  water  from  outboard  flooding  extends  much  farther  than  for  any 
other  type  of  flooding  (Fig.  38).  The  optimum  angle  of  the  outboard  motor  could  be  determined 
very  quickly  during  operation,  and  it  could  be  turned  to  direct  the  discharge  in  different  directions. 
This  method  was  therefore  better  in  providing  uniform  coverage.  The  operator,  however,  was 
required  to  stand  in  ankle  deep  water  during  the  flooding  (Fig.  39)  and  footing  was  treacherous 
unless  crampons  were  used. 

4)  Fire  hoses  were  used  to  distribute  water  from  the  header  into  the  center  of  the  pond.  They 
were  extended  from  the  header  on  the  south  side  of  the  pond  out  onto  the  center  of  the  ice  sheet. 
Three  200-ft  fire  hoses  were  installed  at  60-m  intervals  along  the  header.  The  pressure  loss  in  the 
hoses  was  so  large  that  only  minimal  flow  was  achieved  at  the  outlet.  Erosion  was  also  a big 
problem,  and  an  attendant  was  required  to  change  the  position  of  the  hose  every  half-hour.  The 
wet  2’/i-in.  hose,  partially  frozen  into  the  ice,  was  very  difficult  for  one  person  to  successfully 
wrestle  into  a new  position.  After  5 hours  of  exhausting  effort,  this  method  was  abandoned.  The 
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Figure  39.  Flooding  with  outboard  motor. 


I 


* 


effort  yielded  10  cm  of  flood  water  in  the  center  of  the  sheet;  however,  the  area  covered  was  only 
25%  of  the  pond  surface,  giving  a net  flood  rate  of  2.5  cm  in  5 hours,  or  0.5  cm  hr-1. 

The  above  techniques  were  found  to  be  practical  for  localized  ice  building  only.  They  were  not 
suitable  for  large  scale  coverage  on  a pond  of  this  size.  An  additional  disadvantage  was  that  the 
operator  had  to  remain  out  in  the  cold  for  long  periods,  and  because  of  the  time-consuming  nature 
of  separately  flooding  each  small  area,  the  maximum  depth  that  could  be  obtained  in  any  one  flood 
area  was  7.6  cm  per  8-hour  day  (0.95  cm  hr~! ).  (And  this  was  at  temperatures  of  -20°C,  much 
colder  than  the  average  conditions  considered  in  the  growth  rate  calculations.) 

Ice  volume  measurements 

The  surface  area  of  the  ice  sheet  was  nearly  constant,  varying  only  a few  inches  when  the  pond 
level  changed.  The  volume  was  then  easily  determined  by  measuring  the  thickness.  Thickness, 
however,  varied  widely  due  to  nonuniform  flooding,  and  therefore  a large  number  of  measurements 
were  taken  so  that  averages  could  be  determined  to  give  the  total  volume.  A measurement  net  was 
established  with  lines  at  60-m  intervals  along  the  length  of  the  pond  and  15-m  intervals  across  the 
pond.  This  gave  a total  of  20  positions  for  monitoring  thickness.  Under  some  circumstances  more 
holes  were  drilled  to  determine  local  changes  that  did  not  fit  into  the  average  pattern.  During 
melting,  the  holes  were  redrilled  every  24  hours  to  monitor  the  thickness  loss.  A common  wood 
auger  1 in.  in  diameter,  with  the  center  screw  ground  away  and  an  extension  welded  on  for  greater 
length,  performed  well.  A battery-operated  3/g-in.  electric  drill  was  used  to  drive  the  auger  (Fig.  40). 
Thickness  was  measured  with  a hooked  wire  gage  calibrated  in  inches.  The  wire  was  inserted  in  the 
hole  and  drawn  out  until  the  hook  contacted  the  lower  surface  of  the  ice.  The  thickness  was  then 
read  from  the  scale  on  the  wire  (Fig.  41). 


Melting  experiments 

At  the  onset  of  an  ice  fog  episode,  the  thickness  of  the  existing  ice  sheet  was  monitored  in  the 
above  manner.  The  diversion  of  water  to  the  intake  pond  was  then  stopped  and  the  entire  output 
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of  cooling  water  from  the  power  plant  was  directed  under  the  ice.  The  thickness  was  remeasured 
at  24-hour  intervals,  and  the  amount  of  ice  that  had  melted  in  the  interim  was  determined.  Melting 
of  the  sheet  was  very  dependent  upon  position.  The  first  area  to  melt  was  directly  above  the  input 
header  at  the  west  end  of  the  receiving  pond.  As  soon  as  open  water  was  exposed  above  the  header, 
the  edge  of  the  ice  sheet  melted  from  the  inlet  to  the  distal  end  until  the  entire  pond  was  open  water 
(Fig.  42  and  43).  Melting  on  the  underside  of  the  sheet  showed  a typical  boundary  layer  pattern, 
with  the  heaviest  melt  occurring  in  the  first  60  m from  the  receding  edge  of  the  ice.  Progressively 
less  melting  was  noted  as  the  distance  from  the  edge  increased,  until  at  1 80  m from  the  edge  no 
melting  of  any  consequence  took  place  (Fig.  44).  This  is  to  be  expected  of  flow  along  a quasi- 
infinite  plate  where  a boundary  layer  between  the  moving  fluid  and  the  solid  surface  is  established 
and  continually  thickens  in  the  direction  of  flow.  The  thickening  boundary  layer  provides  progres- 
sively more  and  more  insulation  between  the  warm  water  and  the  ice  sheet.  Heat  transfer  to  the 
sheet,  consequently,  progressively  decreases. 

The  heat  balance  in  the  pond  can  be  expressed  by 

- ~ ~ ~ Os  = *?out 

where  Qin  = heat  input  to  the  pond 
(?E  = evaporative  heat  transfer 
Qc  = conductive  heat  transfer 
(?H  = convective  heat  transfer 
Qr  = radiation  heat  transfer 
<?s  = heat  storage 
(?m  = heat  of  melting 

Oout  = heat  residual  in  water  leaving  the  pond. 

Terms  of  the  equation  have  varying  proportions  depending  on  the  time  since  the  beginning  of  the 
melt  (Table  VIII).  Heat  loss  through  evaporation  <?E  is  initially  zero,  since  the  evaporation  from  the 
ice  surface  is  negligible.  As  open  water  appears,  (?E  starts  to  grow  and  increase  in  direct  proportion 
to  the  open  water  surface  area. 

Total  conductive  heat  loss  Qc  through  the  ice  is  initially  somewhat  greater  than  that  of  evapora- 
tion. It  decreases  as  the  ice  surface  area  decreases  and  is  soon  much  smaller.  The  convective  heat 
loss  Qh  over  the  open  water  is  a function  of  the  air  moving  over  the  pond  and  the  percentage  of 
the  pond  that  is  open  water.  Radiative  heat  loss  from  the  water  (?R  also  appears  as  a small  portion 
of  the  total  heat  loss  as  melting  begins.  It  is  also  a function  of  open  water  and  surface  temperature. 
Radiation  and  convective  losses  over  the  ice  are  contained  in  the  conductive  loss  through  the  ice, 
and  the  assumption  is  made  that  the  radiation  and  convective  losses  are  large  enough  to  remove  all 
of  the  heat  conducted  through  the  ice  sheet.  Storage  C?s  starts  to  increase  as  soon  as  the  warm  water 
is  diverted  back  into  the  receiving  pond.  A component  of  the  storage  term  occurs  due  to  the  influx 
calculations  show  it  to  be  a third-order  term.  Heat  transferred  to  storage  is  large  at  first  but  con- 
tinually decreases. 

The  heat  lost  to  melting  Qm  increases  with  time.  This  is  a function  of  the  temperature  of  the 
ice  sheet,  which  initially  is  low  and  must  be  raised  to  0°C  before  melting  can  proceed.  As  the 
temperature  of  the  ice  sheet  rises,  more  of  the  heat  transferred  to  the  sheet  will  go  into  melting 
and  will  then  show  up  in  the  measurement  of  ice  volume  lost  by  melting. 

Instrumentation  to  measure  the  exact  proportions  of  each  of  these  terms  was  not  available,  nor 
was  the  required  manpower.  They  can,  however,  be  estimated  from  the  measurements  taken  (Table 
VIII).  Heat  lost  by  melting  Qm  was  determined  by  measuring  the  ice  volume  lost  over  consecutive 
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Figure  42.  Melting  progress  in  Eielson  A FB  cooling  pond  ice  sheet,  29-31  January 

and  1-2  February  1973. 
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Figure  44.  Longitudinal  cross  section  of  ice  thickness  during  melting  at  Elelson 
AFB  cooling  pond,  January  and  February  1973. 


Table  VIII.  Heat  loss  components  during  ice  sheet  melting  (W  x 10~s). 


Period 

1 

2 

3 

4 

5 

Average 

Date  (1973) 

1/29-30 

1/30-31 

1/31-2/1 

-2 

2/26-27 

iK-r,) 

20 

33 

38 

28 

23 

Open  water  (m2) 

3000 

3600 

7700 

12600 

3000 

Duration  (hours) 

27 

19 

24.5 

24.5 

24 

<?E 

2.4 

4.6 

11.6 

15.) 

2.5 

7.2 

26.4 

16.2 

17.3 

24.2 

27.7 

22.4 

2.1 

40.2 

5.8 

66.5 

28.6 

3.6 

7.9 

22.3 

25.2 

3.8 

12.6 

Qr 

3.6 

5.8 

14.5 

18.9 

3.7 

9.3 

28.4 

32.2 

74.8 

90.0 

96.1 

64.2 

Mother 

-1.7 

-9.3 

~1.7 

-1.2 

-19.3 

-6.0 

66.1 

54.5 

179.0 

178.0 

181.0 

132.0 

OJQin 

0.42 

0.59 

0.42 

0.51 

0.53 

0.49 

24-hour  periods  from  the  ice  sheet  thickness  measurements  discussed  above.  Evaporation  heat  loss 
Qe  was  determined  from  data  taken  from  the  evaporation  pans  on  the  intake  pond  (see  Section  II). 
Values  for  radiation  heat  loss  from  the  water  surface  QR  are  derived  using  the  considerations 
discussed  in  Section  III.  The  revised  Anderson  equation  was  used  for  atmospheric  back  radiation 
from  the  air  above  the  water. 

Convective  heat  losses  (?H  are  from  values  determined  in  Section  IV.  Storage  Qt  was  determined 
from  measurements  of  the  water  temperature  at  the  circulation  dam  in  the  channel  between  the 
receiving  and  intake  ponds.  Heat  into  the  pond  was  calculated  from  measurements  of  the  flow  rate, 
and  outgoing  temperatures. 
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The  net  flow  rate  into  the  pond  was  measured  by  the  following  method: 

Pond  depth  was  closely  monitored  until  a uniform  depth  was  maintained  for  several  days  so  that 
the  ground  water  recharge  or  discharge  was  negligible.  The  make-up  wells  were  then  activated  until  the 
pond  level  had  risen  5 cm  (2  in.)  above  the  base  level.  At  this  time  both  the  wells  and  the  hot  water 
return  were  stopped.  The  rate  of  fall  of  the  surface  was  timed  until  the  level  had  fallen  to  5 cm 
below  the  base  level.  Input  water  was  then  directed  back  into  the  pond  and  the  level  monitored  to 
check  to  be  sure  that  the  level  returned  to  the  base  line.  The  rate  of  flow  in  and  out  of  the  pond 
could  then  be  determined  with  considerable  accuracy. 

The  heat  input  to  the  pond  (in  Table  VIII)  was  less  during  the  first  two  days  of  melting,  due  to 
diversion  of  a portion  of  *he  input  water  to  a small  creek.  The  average  normal  heat  input  (reflected 
by  the  last  three  terms  of  Table  VIII)  was  1.79x107  W.  A pan  constant  of  1.0  was  used  for  evapora- 
tion (see  Section  II). 

The  surface  areas  of  the  remaining  ice  sheet  and  open  water  could  not  be  measured  with  great 
precision,  since  approaching  the  thin  edge  of  the  melting  ice  was  precarious.  However,  a great  deal 
of  care  was  taken  to  see  that  the  measurements  were  as  accurate  as  possible. 

Ice  fog  suppression  considerations 

Since  nearly  half  of  the  heat  input  to  the  pond  is  used  to  melt  ice  and  most  of  the  rest  is  dissi- 
pated through  radiation,  evaporation,  and  convection,  and  since  it  was  possible  to  maintain  an  ice 
cover  for  several  days,  ice  fog  suppression  with  a manufactured  ice  cover  is  considered  to  be 
feasible.  The  limitations  to  the  technique  are  primarily  those  of  ice  sheet  formation.  If  an  average 
sustained  rate  of  freezing  of  0.75  cm  hr-1  is  maintained,  a 1-m  thick  sheet  of  ice  will  be  formed  in 
5.5  days.  This  volume  of  ice  will  require  1 2 days  to  melt  when  absorbing  50%  of  the  total  heat  in- 
put. A 1 'A-m  sheet  of  ice  would  require  8 days  to  form  and  would  suppress  ice  fog  for  1 8 days. 

(VA  m of  ice  is  the  maximum  thickness  possible  without  considerable  modification  of  the  pond,  and 
1 m is  probably  a more  practical  limit.)  The  longest  period  of  continuous  ice  fog  weather  on  record 
at  the  Fairbanks  Weather  Bureau  is  19  days,  and  66%  of  the  incidents  are  8 days  or  less.  Thus,  a 
1'A-m  sheet  could  suppress  ice  fog  for  even  the  longest  periods,  and  a 1-m  thickness  would  suffice 
for  the  shorter  ones. 

Between  ice  fog  sessions,  there  must  be  sufficient  time  to  rebuild  the  ice  cover.  The  situation 
of  a long  ice  fog  experience,  followed  by  a second  similar  period  with  only  1 or  2 days  between, 
would  require  the  thicker  1!4-m  sheet.  The  average  time  interval  between  ice  fog  incidents  in  the 
last  10  years  is  6.2  days.  This  would  allow  adequate  time  to  rebuild  a 1-m  thick  ice  sheet  at  a rate 
of  0.75  cm  hr'1 . On  26  occasions  during  the  last  10  winters,  a break  of  2 days  or  less  fell  between 
ice  fog  incidents,  but  on  only  two  of  these  occasions  was  the  total  sequence  of  ice  fog  events  long 
enough  to  exhaust  a 1-m  thick  ice  sheet. 

During  the  melting  of  the  sheet,  ice  fog  will  be  generated  as  soon  as  the  first  open  water  appears; 
however,  the  percentage  of  the  normal  production  of  fog  is  very  small.  If  suppression  is  defined  as 

5 = evaporation  during  suppression 
evaporation  without  suppression 

the  graph  (Fig.  45)  shows  the  performance  during  the  tests  of  this  experiment  in  January  and 
February  1973.  Figures  42  and  43  show  the  limits  of  the  ice  sheet  and  the  ice  thicknesses  measured 
during  two  melting  experiments.  Figure  46  and  the  cover  photograph  show  the  pond  with  and  without 
complete  ice  cover.  Figures  47-50  show  events  and  problems  of  relevance  during  the  experiments. 
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Figure  45.  Ice  fog  produced  during  suppression,  January  and  February  1973. 


Figure  46.  Ice  fog  suppression  with  total  ice  cover  (~30°C). 


Experimental  results 

The  unusually  mild  winter  of  1973,  coupled  with  the  inability  of  the  installation  to  form  a 
uniformly  thick  ice  sheet,  allowed  only  two  melting  experiments.  In  addition,  an  ice  sheet  thick- 
ness of  53  cm  was  the  maximum  obtained.  This  is  only  half  of  the  thickness  felt  to  be  optimum. 
The  maximum  long-term  freezing  rate  obtained  during  the  ice  building  experiments  was  7.5  cm  in 
24  hours  {0.31  cm  hr-1 ).  The  basic  limitation  was  that  of  freezing  from  the  top  down.  No  way 
was  available  to  prevent  the  water  from  puddling  and  forming  a skim  of  ice  over  the  top.  As  soon 
as  this  occurred,  all  further  freezing  was  then  done  by  transferring  the  heat  through  the  ice  on  top 
of  the  puddles,  in  accordance  with  the  limitations  of  eq  60.  This  was  the  fundamental  freeze  rate 
limitation. 
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Figure  50.  Ice  breakup  during  final  stages  of  melting. 
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The  optimum  practical  thickness  of  each  applied  layer  of  water  to  yield  a maximum  growth  rate 
is  a function  of  the  temperature  and  weather  conditions,  as  well  as  the  method  of  application  used. 
An  overall  average  of  0.46  cm  hr-1  of  ice  growth  was  measured  for  short-term  experiments  during 
the  winter.  Since  the  growth  rate  is  proportional  to  the  square  root  of  time,  the  highest  rate  of 
freezing  occurs  during  the  first  few  minutes  after  the  film  of  ice  is  formed.  Thin  films  of  water 
would  therefore  freeze  much  faster  than  the  0.46  cm  film.  If  a very  thin  film  could  be  applied  at 
frequent  intervals,  the  overall  result  might  well  yield  an  optimum  growth  rate  of  perhaps  0.75  cm 
hr'1 , approaching  the  theoretical  maximums  discussed  earlier. 

Although  the  above  technique  appears  to  offer  a workable  means  of  ice  fog  suppression,  the 
limitation  of  ice  growth  rate  may  well  be  an  overriding  consideration  that  limits  its  practicality. 


SECTION  VII.  RECOMMENDATIONS  AND  CONCLUSIONS 

The  work  described  in  this  report  has  tested  the  feasibility  of  two  techniques  of  ice  fog  suppres- 
sion. Conclusive  tests  have  not  been  made  to  determine  which  is  the  most  practical  or  efficient. 
Further  studies  are  needed,  beginning  with  more  exhaustive  tests  of  hexadecanol.  In  view  of  the 
present  work,  it  appears  that  hexadecanol  may  provide  the  most  inexpensive  and  efficient  solution 
of  the  problem.  Should  further  hexadecanol  experiments  show  this  technique  to  be  impractical,  a 
decision  would  have  to  be  made  concerning  the  desirability  of  continued  use  of  the  pond  during 
the  winter  months,  as  opposed  to  the  injection  well  approach  discussed  in  Section  V.  If  an  impact 
study  of  the  injection  well  technique  shows  that  there  will  be  no  detrimental  effects  to  the  operation 
of  the  power  plant,  and  the  capital  expenditures  are  within  reason,  then  this  is  the  next  most  logical 
approach.  If  the  implementation  of  an  injection  well  proves  impractical,  the  installation  for  latent 
heat  storage  suppression  appears  to  offer  the  most  promise. 

Hexadecanol  studies 

The  fundamental  questions  to  be  resolved  concerning  this  approach  are  the  following: 

1 ) Film  integrity.  A floating  grid  or  net  should  be  fabricated  and  tested  for  effectiveness.  The 
grid  size  and  height  of  protrusion  above  water  need  *o  be  optimized.  This  will  provide  a means  for 
the  film  to  maintain  continuity  over  the  entire  pond  and  eliminate  the  problem  of  wind  concentrating 
the  film  on  one  side  of  the  pond. 

2)  Film  life.  Biological  degradation  is  of  concern  in  relation  to  film  life.  The  extent  of  the 
problem  should  be  determined.  Methods  of  eliminating  or  reducing  its  magnitude  by  poisoning  the 
film  with  a bactericide  or  treating  the  pond  between  applications  of  the  alcohol  need  to  be  studied 
and  techniques  developed.  Film  life  should  be  determined  under  ice  fog  conditions. 

3)  Alcohol  blends.  Several  other  long-chain  fatty  alcohols  need  to  be  tested  to  determine  which 
one,  or  which  blend  of  two  or  more,  gives  the  best  suppression  of  evaporation  and  longest  lifetime. 

4)  Application  technique.  Many  methods  have  been  used  to  apply  the  alcohols  to  the  water  sur- 
face. Application  can  be  done  with  a hand-held  lawn  spreader  as  was  done  in  the  studies  of  this 
report,  or  by  spraying  slurries  of  alcohol  and  water.  A method  compatible  with  the  conditions  of 
ice  fog  needs  to  be  developed  and  implemented. 

The  tests  thus  far  are  not  conclusive,  but  strongly  suggest  that  alcohol  application  is  an  effective 
suppression  technique.  Suppression  of  84%  at  temperatures  below  freezing  appears  to  be  better  than 
the  62%  reported  under  warmer  conditions.  The  major  problem  will  probably  be  that  of  maintaining 
film  integrity  during  wind.  If  practicality  can  be  demonstrated,  use  of  hexadecanol  may  be  the  least 
expensive  method  by  a wide  margin.  Application  should  pose  no  technical  problems. 


Injection  well  suppression 

Before  this  method  can  be  used,  several  questions  need  to  be  resolved.  They  are  the  following: 

1 ) Cooling  water  supply.  Tests  must  be  made  to  determine  if  the  make-up  wells  can  supply  ade- 
quate volumes  of  water  on  a continuous  basis.  Tests  run  during  the  winter  indicate  that  the  wells 
are  only  producing  at  approximately  half  of  their  original  capacity. 

2)  Raw  water  problems.  The  detrimental  effects  of  using  raw  water  in  the  condensers  need  to  be 
studied  and  any  problems  resolved. 

3)  Injection  well.  Injection  wells  are  known  for  gradually  building  up  resistance  to  water  injec- 
tion due  to  "plugging.”  A well  or  wells  of  adequate  size  must  be  developed  and  tested  along  with 
an  emergency  back-up  system.  (A  valving  arrangement  to  allow  the  cooling  water  to  be  changed 
back  to  the  cooling  pond  would  provide  adequate  safety.) 

4)  Recirculation  bypass.  A method  of  recirculating  a portion  of  the  warm  cooling  water  needs 
to  be  designed  into  the  system  to  allow  the  cold  well  water  to  be  tempered  before  it  is  sent  into 
the  condensers. 

5)  Costs.  The  initial  cost  of  the  well  drilling,  injection  pumps,  internal  piping,  and  valving  needs 
to  be  determined  to  assure  that  this  method  is  economically  feasible. 

The  present  well  installation  at  Eielson  is  inadequate  for  direct  use  of  this  method  but  it  could 
be  sufficiently  enlarged.  The  injection  of  heated  water  into  wells  provides  a positive  elimination  of 
the  ice  fog  contribution  of  the  ponds;  however,  it  poses  several  environmental  questions. 

Latent  heat  storage  suppression 

The  critical  limiting  factor  in  this  technique  is  the  ability  to  rapidly  form  a uniform  ice  sheet  of 
sufficient  thickness.  One  means  is  to  install  spray  nozzles  in  a 30-m  grid  so  that  the  ice  sheet  can 
be  uniformly  covered  by  spray.  This  could  be  accomplished  with  three  headers  equally  spaced 
across  the  width  and  running  the  length  of  the  pond.  Nozzles  capable  of  spraying  1 5 m throughout 
360°  would  be  installed  every  30  m.  The  headers  could  be  installed  either  under  the  surface  of  the 
pond  (and  be  provided  with  an  air  purge  so  that  the  nozzle  standpipe  would  not  freeze  in  the  area 
where  it  passes  through  the  ice  sheet),  or  alternatively,  they  could  be  set  on  pilings  over  the  pond 
(with  the  nozzles  on  the  bottom  of  the  pipe).  During  the  experiments  discussed  in  Section  VI,  no 
freezing  problems  were  experienced  after  the  larger  (1 ,9-cm-diameter)  nozzles  were  installed. 

Subject  to  this  consideration,  the  technique  has  proven  to  be  a practical  working  means  of  ice 
fog  suppression.  Visibility  conditions  in  the  vicinity  of  the  pond  were  dramatically  improved  during 
ice  fog  weather  when  suppression  was  in  effect.  No  discernible  fog  could  be  seen  during  this  period, 
and  bushes  were  visible  in  excellent  detail  across  the  length  of  the  pond.  Some  problems  need  to 
be  considered,  however,  in  the  application  of  this  technique: 

1.  A circulation  dam  between  the  two  sides  of  the  pond  must  be  installed  to  eliminate  the  density 
gradient  circulation  from  the  warm  side  to  the  cold  side  of  the  pond.  This  provides  the  isolation 
needed  to  form  an  ice  sheet  at  the  maximum  rate  in  all  weather  conditions. 

2.  To  provide  better  control  of  the  temperature  of  the  intake  pond,  the  existing  make-up  wells 
should  be  rerouted  to  empty  into  either  side  of  the  pond.  This  will  enhance  operating  efficiency 
and  reduce  time  required  to  switch  from  an  ice  building  to  an  ice  suppression  configuration. 

Comparison  of  approaches 

Of  the  three  above  recommended  approaches,  the  hexadecanol  appears  to  have  the  advantages 
of  low  initial  cost  and  low  maintenance  and  operating  cost.  The  injection  well,  while  holding  a 
good  deal  of  promise,  needs  careful  study  before  the  rather  formidable  initial  investment  is  justified. 
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The  latent  heat  storage  technique,  while  appearing  to  be  the  third  choice,  has  the  big  advantage  of 
having  been  proven  in  the  field  under  ice  fog  conditions.  We  know  more  about  this  method  and  its 
requirements  and  performance  than  any  of  the  others  at  this  time. 

Measurements  of  evaporation 

Evaporation  into  below-freezing  air  averages  5 kg  m-2  day"' . Evaporation  measurements  are 
difficult  under  the  very  cold  conditions  of  ice  fog,  and  a scatter  of  data  points  is  an  inevitable  re- 
sult. However,  the  data  were  of  sufficient  quantity,  and  scatter  of  the  same  magnitude  as  the  con- 
trolled laboratory  experiments  of  Yen  and  Landvatter  (1970)  was  displayed.  This  lends  confidence 
to  the  average  values  obtained  in  these  studies.  The  formulas  developed  from  this  work  should  be 
useful  for  ice  fog  evaporation  calculations,  and  their  accuracy  should  be  an  improvement  over  those 
equations  previously  available. 

Measurements  of  radiation  loss 

Longwave  radiation  losses  are  a significant  portion  of  the  total  heat  losses.  They  are  attenuated 
by  the  ice  fog  cover,  but  even  with  heavy  ice  fog,  they  continue  to  be  significant.  Average  net 
longwave  radiation  lost  from  the  surface  was  1 1 9 W m~2. 

Convective  heat  losses 

Convective  heat  losses  were  estimated  using  both  Rimsha  and  Donchenko’s  equation  (1958), 
Kays’  equation  (1966),  and  Rotem  and  Claassen’s  equation  (1969).  Agreement  between  the  various 
equations  was  very  good  at  lower  temperature  differences,  but  during  periods  with  low  winds,  they 
diverged  considerably  as  temperature  differences  increased.  More  work  is  needed  to  verify  the 
correct  convective  heat  loss  equations  at  these  times. 

Ice  fog  suppression 

The  prime  purpose  of  this  paper  was  the  investigation  of  techniques  for  the  suppression  of  ice 
fog  from  open  water  of  cooling  ponds.  Values  for  the  various  modes  of  heat  transfer  were  measured 
during  the  studies.  These  values  were  used  for  engineering  estimates  of  the  neat  budget  during  the 
melting  phase  of  the  experiments  on  latent  heat  suppression.  Equations  were  developed  to  make  it 
possible  to  determine  the  necessary  heat  balance  components  of  the  pond. 

In  the  past,  ice  fog  has  been  tolerated,  studied,  and  cursed.  No  real  attempt  has  been  made  to 
suppress  it  at  its  source.  It  is  the  conclusion  of  this  paper  that  suppression  is  not  only  feasible  - but 
feasible  by  several  techniques.  Suppression  both  by  latent  heat  storage  and  by  monomolecular  films 
was  proven  to  be  effective.  In  an  engineering  assessment  of  the  problem,  both  the  effectiveness  and 
the  economics  of  each  method  must  be  considered.  From  a cosf  standpoint,  the  monomolecular 
film  technique  appears  to  be  better.  From  the  standpoint  of  effectiveness,  the  latent  heat  storage 
method  was  somewhat  more  effective  and  more  positive.  Either  way,  ice  fog  from  open  water  sur- 
faces can  be  suppressed  by  using  the  engineering  techniques  discussed.  The  reduction  in  the  total 
ice  fog  problem  should  be  significant,  and  will  be  visually  evident  to  those  working  and  living  with 
the  problem. 
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APPENDIX  A.  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT  FROM 
WATER  TO  ICE  AND  FROM  ICE  TO  AIR 

Kays  (1966)  develops  the  solution  for  turbulent  flow  of  a fluid  past  a semi-infinite  plate  of 
constant  temperature.  This  meets  the  conditions  of  the  cooling  pond  at  Eielson  Air  Force  Base  if 
we  restrict  the  solution  to  low  velocity  flow,  constant  property  fluids,  and  Prandtl  numbers  that 
range  from  0.5  to  10.  These  restrictions  are  easily  met  for  the  case  in  question  and  the  general 
formula  developed  should  give  a good  comparison  for  measured  values. 

Turbulent  heat  transfer  in  the  Prandtl  number  Pr  range  of  interest  can  be  expressed  by  the  rela- 
tion: 

Stx  Pro.4  = 0.0295  Re-0-2  (A1) 

where  the  Stanton  number  St  is 


St. 


PC p 


(A2) 


and  x is  the  distance  from  the  leading  edge  of  the  plate.  Rearranging  and  substituting  for  Stx  yields 
the  convective  heat  transfer  coefficient  hx.  Then 


h = (0.0295)(p  Cp  Q 
x (Pr°-4)(Re0-2) 

where  Cp  = specific  heat 
p = density 
u = velocity  of  water 
Re  = Reynolds  number. 

Evaluating  this  relation  for  water  at  the  following  conditions: 


T„  = 10°C 

u„  = 0.152  (m  min'1) 

where  T„  = free  stream  temperature 
um  = free  stream  velocity 

gives 

Pr  = 7.0 
Cp  = 1.0 

v - 7.8  x10-5  (m2  s'1) 
p = 999  (kg  m'J) 
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where  v = the  kinematic  velocity. 


r.  > = 174,9,9 

" 7.8xl0*s 

where  x = distance  along  ice  sheet. 

Then 


(0.0295)(pCo</J(1000) 

— — — — 1-  = 1276  Wm-2°C-‘ 

(7)0-4  (174.919)0-2 


Evaluating  the  relation  for  air  at  the  following  conditions: 
r„  = -20  (°C) 

= 1-0  (m  s'1) 
x = 1 (m) 

I gives 

F 

\ 

[ Pr  - 0.72 

; Cp  = 0.24  (cal  g'>  0C'«) 

| v = 0.93x10“s  (m2  s'1) 

: Re  = 107,527. 

■ Then: 

h = (0-0295)  (1 .4  x10r3)(0.24)(l  x106)(1.0)(4.184) 
*air  (0.72)0-4  (T  .08  x TO5)0-2 

t 

F 

i /»,  = 4.66  Wm'2  °C-'. 
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